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CHAPTER   I 
INTRODUCTION 

The  development  of  design  and  analysis  techniques  for 
plate  cams  has  remained  an  active  research  area  for  several 
decades.  There  are  many  books  devoted  entirely  to  cam 
design  [1,  2],  while  most  machine  design  books  have  at  least 
one  chapter  on  this  subject.  Numerous  papers  on  this  topic 
are  still  published  every  year,  which  indicates  that  a  lot 
of  research  still  remains  to  be  done  in  this  area. 

The  use  of  standard  curves  is  the  traditional  approach 
to  cam  design.  The  most  commonly  used  curves  are  harmonic 
and  cycloidal  curves.  In  this  approach,  the  design  process 
is  simply  the  selection  of  a  set  of  curves  which  satisfies 
the  required  boundary  conditions.  The  process  of  selecting 
the  curves  is  a  tedious  procedure,  since  different  sets  of 
curves  are  required  for  different  design  cases.  Moreover, 
this  procedure  restricts  the  range  of  solutions,  since 
standard  curves  are  only  a  small  subset  of  the  possible 
solutions  which  satisfy  the  specified  boundary  conditions. 

Another  common  approach  in  designing  cams  is  to  use 
polynomial  curves  to  represent  the  cam  motions  [ 3 ] .  Depend- 
ing on  the  boundary  conditions,  a  polynomial  curve  with  the 
number  of  coefficients  equal  to  the  number  of  boundary 
conditions  is  generated.   Therefore,  the  solution  to  a 
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motion  specification  is  uniquely  determined.  Obviously, 
both  standard  curves  and  polynomial  curves  only  permit  the 
design  to  satisfy  the  boundary  conditions;  thus  the  solution 
may  not  satisfy  other  criteria  related  to  performance  and 
manufacture. 

An  improved  approach  is  presented  in  [4].  In  this  case 
also,  polynomial  curves  are  used  to  design  the  cam.  Howev- 
er, the  order  of  the  polynomial  curve  is  set  to  one  order 
higher  than  that  required  to  satisfy  the  boundary  condi- 
tions. By  doing  so,  a  free  coefficient  is  generated,  and  it 
is  treated  as  a  design  variable.  This  design  variable  can 
be  used  to  satisfy  performance  and  manufacturing  aspects 
through  an  optimization  process.  More  free  coefficients  can 
be  generated  by  raising  the  order  of  the  polynomial  curve. 
However,  there  is  a  practical  limit  on  this  approach,  since 
a  high  order  polynomial  curve  causes  undesirable  oscilla- 
tions in  the  cam  motion.  Besides  the  limited  number  of  free 
coefficients,  polynomial  curves  do  not  provide  good  local 
control  on  cam  motion  and  thus  may  significantly  limit  the 
solution  to  a  design  problem.  Furthermore,  due  to  the  lack 
of  local  control,  this  technique  may  not  be  able  to  provide 
any  solutions  for  problems  where  rigorous  performance  re- 
quirements are  to  be  enforced  in  the  interior  of  a  rise  or 
return  segment. 


The  purpose  of  this  work  is  to  present  a  new  design 
approach  which  is  more  powerful  than  currently  available 
design  methods.  The  approach  adopted  is  to  use  B-spline 
curves  to  design  the  cam.  An  example  of  a  B-spline  curve  is 
shown  in  Figure  1.1. 

Some  preliminary  work  on  the  application  of  B-splines 
in  cam  design  has  been  done  [5].  However,  the  overall 
development  in  this  area  is  very  limited,  and  automated 
design  procedures  are  practically  non-existent. 

B-spline  curves  are  parametric  piecewise  polynomials 
which  are  defined  by  a  series  of  control  points.  A  B-spline 
curve  provides  good  local  control  over  its  shape.  Therefore, 
it  can  be  used  to  model  a  large  variety  of  complex  curves, 
which  cannot  be  modeled  by  either  the  standard  curves  or  the 
polynomials.  Since  B-splines  are  conveniently  and  complete- 
ly defined  by  their  control  points,  the  design  of  a  B-spline 
cam  is  a  process  of  placing  the  control  points  at  the  appro- 
priate locations  to  meet  the  specified  requirements. 

In  the  current  industrial  design  environment,  design 
tasks  are  usually  performed  without  considering  the  relevant 
manufacturing  aspects.  Without  proper  communication  between 
design  and  manufacture,  the  solution  obtained  by  this  ap- 
proach may  not  be  one  that  can  be  easily  produced.  Efforts 
have  been  made  to  link  the  design  and  manufacturing  stages 


a  Control  point 


Figure  1.1  :  A  B-spline  Curve 


more  tightly.  These  considerations  provide  a  strong  motiva- 
tion to  consider  both  the  design  and  manufacturing  aspects 
simultaneously,  so  that  the  final  design  will  satisfy  the 
requirements  of  both  aspects. 

In  this  thesis,  a  computer-oriented  method  for  inte- 
grating the  design  and  manufacture  of  B-spline  plate  cams  is 
developed.  This  approach  allows  the  designer  to  enforce  not 
only  the  performance  requirements,  but  also  the  manufactur- 
ing concerns.  In  this  environment,  the  designer  may  impose 
various  performance  requirements  on  quantities  such  as  the 
pressure  angle,  contact  stress,  follower  acceleration,  etc., 
as  well  as  on  manufacturing  considerations,  such  as  the 
cutting  time,  material,  etc.  This  results  in  the  establish- 
ment of  a  stronger  link  between  design  and  manufacturing. 

In  the  present  work,  the  linking  of  design  and  manufac- 
turing aspects  is  achieved  by  using  nonlinear  optimization 
techniques.  In  order  to  use  the  techniques  of  optimal 
design  to  solve  a  given  problem,  the  problem  must  be  stated 
in  a  standard  format.  The  format  of  the  standard  nonlinear 
programming  problem  (NLP)  is  the  following: 

Suppose  the  design  of  a  system  is  specified  by  a  vector 
b  of  r  design  variables: 

b  =  [bl,b2,b3, ,br]T 


Find  the  value  of  b  to  minimize  a  given  cost  function,  F(b) , 
subjects  to  the  constraints 

gj(b)  <  0  j  =  l,m 

b-l(b)  =0  1-1,11 

The  advantage  of  implementing  NLP  is  that  it  allows 
integration  of  design  and  manufacture.  Based  on  the  speci- 
fied requirements,  integration  can  be  achieved  by  selecting 
the  design  variables  from  both  aspects.  Similarly,  the 
constraint  functions  can  also  be  derived  from  both  design 
and  manufacturing.  Once  the  optimization  problem  is  formu- 
lated, there  are  many  available  algorithms  for  solving  the 
problem  automatically  [6,  7].  As  a  result,  implementing  NLP 
in  cam  designs  not  only  reduces  the  amount  of  designer  time 
needed,  but  also  reduces  the  level  of  skill  required  to 
produce  the  design. 

Two  different  aspects  of  manufacturing  are  considered 
in  this  thesis:  1)  the  cutting  time  to  manufacture  the  cam, 
and  2)  the  material  required  to  produce  the  cam.  In  the 
first  case,  it  is  assumed  that  the  cam  will  be  manufactured 
on  an  NC  machine  by  an  end  milling  process.  It  is  also 
assumed  that  the  cam  will  be  cut  from  a  rectangular  blank, 
once  the  dimensions  of  the  cam  are  obtained  from  the  optimi- 
zation procedure.  The  goal  is  to  design  a  cam  which  re- 
quires minimum  machining  time,  and  satisfies  the  specified 
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design  requirements.  Therefore,  the  cutting  time  is  consid- 
ered as  the  cost  function  in  the  NLP.  The  final  design 
obtained  from  the  optimization  process  is  then  utilized  to 
generate  the  NC  codes  to  produce  the  cam  on  an  NC  machine. 

For  the  second  case,  we  consider  the  design  of  a  mini- 
mum size  cam.  In  this  case,  the  cross-sectional  area  of  the 
cam  is  treated  as  the  cost  function. 

The  idea  of  using  optimization  techniques  in  the 
design  of  cams  is  also  applicable  for  problems  where  manu- 
facturing considerations  are  not  the  main  concern.  The 
same  approach  can  be  extended  to  problems  where  the  objec- 
tive function  and  constraints  are  all  drawn  from  design 
considerations. 

Note  that  since  B-spline  curves  are  entirely  defined  by 
their  control  points,  the  coordinates  of  the  points  can  be 
treated  as  design  variables  in  an  optimization  process  for 
the  satisfaction  of  design  requirements.  Further  discus- 
sions on  this  point  are  presented  in  Chapter  II  and  V. 

A  discussion  on  the  design  of  cams  using  B-spline  curves 
is  presented  in  Chapter  II.   Chapter  III  is  devoted  entirely 
to  a  discussion  of  the  design  aspects,  while  Chapter  IV 
presents  the  manufacturing  aspects  related  to  plate  cams. 
In  Chapter  V,  the  formulation  of  the  nonlinear  optimization 


problem  is  presented,  with  emphasis  on  the  selection  of 
design  variables,  cost  functions  and  constraint  functions. 
Chapter  VI  deals  with  the  design  sensitivity  analysis, 
while  Chapter  VII  presents  some  numerical  examples  that  were 
solved  to  verify  the  proposed  technique  and  its  implementa- 
tion. Conclusions  drawn  from  this  work  and  recommendations 
for  future  work  in  this  area  are  summarized  in  Chapter  VIII. 


Chapter  II 
B-SPLINE  CAM  DESIGN 

A  variety  of  mechanical  systems  are  available  for  use 
by  the  designer  for  the  satisfaction  of  specified  motion 
requirements.  The  designer  must  generally  choose  between 
constrained  linkages,  cam  systems  and  control  systems.  Cam 
systems  are  capable  of  generating  a  large  variety  of  output 
motions.  The  flexibility  offered  during  the  design  process 
is  one  of  the  key  advantages  of  cam  systems  [3].  Cams  can 
assume  a  wide  variety  of  physical  shapes,  and  may  drive  many 
different  types  of  followers. 

2.1  Types  of  Cams  and  Followers 

Cams  are  classified  in  terms  of  their  physical  shape. 
Figure  2.1  illustrates  three  different  types  of  cams  : 

a)  A  plate  cam,  also  called  a  radial  cam  or  a  disk  cam. 

b)  A  wedge  cam. 

c)  A  cylindrical  cam  or  barrel  cam. 

Followers  in  general  may  be  classified  based  on   one  or 
more  of  the  following  criteria: 

1)   The  construction  of  the  surface  in  contact,  e.g.,   a 
knife  edge,  roller,  or  flat-face. 
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2)  The  type   of  movement,   whether  translating   or 
oscillating. 

3)  The  location  of  the  line  of  movement  with  reference 
to  the  cam  center,  whether  radial  or  offset. 

Figure  2.2  shows  plate  cams  with  four  different  types  of 
followers. 

This  research  focuses  on  plate  cam  design  with  offset 
translating  roller  follower.  Figure  2.3  shows  a  cam-follow- 
er system  of  this  type.  In  this  cam  system,  the  follower 
translates  along  its  axis,  which  is  displaced  from  the  cam 
center  of  rotation. 


2.2   Lift  Diagram 

A  cam  system  is  usually  driven  by  an  input  shaft  at  a 
constant  velocity  to  generate  the  output  motion  of  the 
follower.  The  input  and  output  motion  of  a  cam  system  can 
be  shown  on  a  lift  diagram.  In  Figure  2.4,  the  abscissa 
represents  one  cycle  of  input  motion.  The  ordinate 
represents  the  output  motion  Y.  The  output  motion  can  be 
related  to  the  input  motion  in  functional  form  : 


X  =  f(S)  (2.1] 
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Figure  2.3  :  Plate  Cam  with  Offset  Translating 
Roller  Follower 
(  Extracted  from  Reference  [3]  ) 


13 


8 

u 

tr- 
io 


0) 

M 

3 

5 


a 
o 
u 

13 

0) 

u 

« 

M 
-u 
X 
H 


where 

Y  is  the  follower  position 

9  is  the  cam  angle  of  rotation 

The  cam  motion  can  be  classified  into  three  regions. 
The  region  where  the  follower  is  moving  away  from  the  cam 
center  is  called  a  rise.  The  maximum  rise  is  called  the 
lift.  A  return  is  a  region  where  the  follower  is  moving 
towards  the  center  of  cam.  A  dwell  is  a  region  in  which  the 
follower  is  at  rest.  In  this  thesis,  it  is  assumed  that  the 
cam  motion  consists  of  single  rise  and  single  return  with 
top  and  bottom  dwells. 

Generally,  the  follower  position  is  normalized  by  divid- 
ing the  actual  position,  Y,  by  the  total  lift.  The  range  of 
the  normalized  follower  position,  denoted  by  y,  is  thus  from 
0  to  1.  The  actual  position  can  be  obtained  by  multiplying 
the  normalized  position  by  the  total  lift  : 


Y  =  L  *  y 

where 

Y  is  the  actual  position 

y  is  the  normalized  position 
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(2.2) 


Likewise,  the  derivatives  of  the  follower  position  with 
respect  to  the  angle  of  rotation  "9"  are  as  follows  : 

Y'  =  L  *  dy/d©  (2.3) 

y"  =  l  *  d2y/d2e  (2.4) 

Y' ' '  =  L  *  d3y/d39  (2.5) 

The  quantities  on  the  left  hand  side  of  Equations  (2.3), 
(2.4)  and  (2.5)  are  also  known  as  reduced  velocity,  reduced 
acceleration  and  reduced  jerk  respectively.  The  true  veloc- 
ity, acceleration  and  jerk  can  be  obtained  by  the  following 
equations  : 


(2.6) 
(2.7) 
(2.8) 


Y  =  Y'  *  u 

Y  =  Y  •  '  *  w  2 
Y*  =  Y'  '  '  *  w3 

where  u  is  the  angular  velocity  of  cam. 


2.3   Boundary  Conditions 

In  the  last  section,  it  was  noted  that  a  complete  cycle 
of  the   cam  motion  consists  of  different  zones.   When  making 
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the  transition  from  one  zone  to  the  next,  the  designer  must 
ensure  continuity  of  position  and  its  derivatives  up  to  the 
degree  required.  This  is  usually  done  by  specifying  appro- 
priate boundary  conditions  at  both  ends  of  the  rise  and 
return  zones.  These  boundary  conditions  are  usually  speci- 
fied in  terms  of  normalized  position,  reduced  velocity, 
reduced  acceleration,  reduced  jerk,  etc.  Figure  2.5  illus- 
trates the  following  boundary  conditions  for  a  rise  zone  and 
a  return  zone: 

Rise   Zone 

8/jS   =    0 

Y  =    0,  Y'    =    0,  Y"    =    0,       Y"  •    =    0 
9//S  =   1 

Y  =    1,  Y'    =    0,  Y 


I    !    1 


Return  Zone 

e//3  =  0 

Y  =   1,  Y'    =   0,         Y'"    =   0 

9//3   =   1 

Y  =    0,  Y'    =    0,  Y"    =    0,  Y1  •  •    =    0 

where  0  is  the  range  of  cam  rotation  for  the  respective 
zone  . 
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Position     Y 

V«Hootty      Y' 

Aeflalar      Y" 
Jark               Y"' 

/   Y"/     V/           >\                     \ 

RISE      CYC 

:i_e 

Figure  2.5(a):   Illustration  of  Boundary  Conditions 
for  A  Rise  Cycle 
(  Extracted  front  Reference  [4]  ) 
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Position     Y 

V.loolty       Y" 

Jark               Y"1 

/                     Xy111 

RETURN 

CYCLE 

Figure  2.5(b):   Illustration  of  Boundary  Conditions 
for  A  Return  Cycle 
(  Extracted  from  Reference  [4]  ) 
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2.4   Techniques  of  Design 

A  variety  of  graphical  and  analytical  techniques  are 
available  to  design  a  cam  system.  In  conventional  design, 
basic  curves  are  used  to  design  the  cam  profile  for  speci- 
fied follower  motions  and  boundary  conditions.  Some  of  the 
basic  curves  that  are  commonly  used  are  cycloidal  curves, 
simple  harmonic  curves,  modified  harmonic  curves, parabolic 
curves,  etc.  Based  on  the  boundary  conditions,  the  designer 
must  select  a  set  of  appropriate  curves  to  meet  these  condi- 
tions. It  is  obvious  that  different  types  of  curves  are 
required  for  different  cases.  Therefore,  designing  a  cam 
system  using  basic  curves  is  a  tedious  procedure.  Further- 
more, this  technique  is  purely  design  for  satisfying  the 
boundary  conditions.  It  does  not  provide  any  control  over 
the  motion  in  the  interior  of  a  region.  This  may  result  in 
serious  problems  such  as  excessive  pressure  angle  and 
stress. 

The  range  of  possible  solutions  can  be  increased  by 
using  polynomials  to  design  the  cam  profile.  This  approach 
can  conveniently  generate  a  wide  variety  of  cam  profiles, 
which  cannot  be  modeled  using  basic  curves.  To  satisfy  the 
boundary  conditions,  a  suitable  order  of  polynomial  is 
required.  The  designer  can  generate  free  variables  by 
raising  the  order  of  the  polynomial  to  an  order  higher  than 
that  required  to  satisfy  the  boundary  conditions  [4].   The 
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Figure  2.6(a):   Displacement  Diagram  and  Derivatives  for 
Full-Rise  Simple  Harmonic  Motion 
(  Extracted  from  Reference  [3]  ) 


Figure  2.6(b):   Displacement  Diagram  and  Derivatives  for 
Full-Rise  Cycloidal  Motion 
(  Extracted  from  Reference  [3]  ) 
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Figure  2.6(c):   Displacement  Diagram  and  Derivatives  for 
Full-Rise  Modified  Harmonic  Motion 
(  Extracted  from  Reference  [3]  ) 


Figure  2.6(d):   Displacement  Diagram  and  Derivatives  for 
Full-Return  Modified  Harmonic  Motion 
(  Extracted  from  Reference  [3]  ) 
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Figure  2.6(e):   Displacement  Diagram  and  Derivatives  for 
Full-Return  Simple  Harmonic  Motion 
(  Extracted  from  Reference  [3]  ) 


Figure  2.6(f):   Displacement  Diagram  and  Derivatives  for 
Full-Return  Cycloidal  Motion 
(  Extracted  from  Reference  [3]  ) 
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extra  coefficients  can  be  used  to  satisfy  more  rigorous 
constraints,  and  to  reduce  cost.  However,  a  very  high  order 
polynomial  will  result  in  oscillation  of  the  cam  profile. 
Therefore,  the  designer  is  restricted  to  a  small  number  of 
free  variables.  Furthermore,  due  to  oscillation,  this 
method  may  be  impractical  for  design  problems  with  several 
boundary  conditions. 

A  different   approach   is   adopted   in  this   research.      Here, 
B-spline  curves   are     used  to  design  the  cam. 

2.5     B-spline  Curves 

B-spline   curves   are   parametric   piecewise   polynomials 
which  are  defined  as   follows    [8]    : 


n 
P(u)    =  S  Pi  Ni  k(u)  (2.9) 

i=0 


where 


u  is  the  independent  parameter 

Pi  is  the  position  vector  of  the  itn  control  polygon 

vertex 
Ni,k(u)  are  blending  functions 
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P(u)   is  the  position  vector  of  the  point   on  the   B- 
spiine  curve  corresponding  to  a  given  value  of  u 

For  a  2-dimensional  problem,  P(u)  =  [  Px(u),  Py(u)  ], 
and  Pi  =  [  Pix,  Piy  ]. 

The  blending  functions,  Ni(lc(u)  are  defined  recursively 
by  the  following  equations  : 

Ni,l(u)  =1   if  ti  <  u  <  ti+1 

(2.10) 
=  0   otherwise 

and 

<u  -  tp  Ni)k_1(U)  (ti+k  -  u)  m+i,k-i(u) 

Ni,k(u)= +  ■ - 

ti+k-1  "  ti  ti+k  -  ti+1 

(2.11) 

The  ti  are  called  the  knots  of  the  B-spline.  They 
represent  the  values  of  the  independent  variable,  u  at  which 
the  curve  switches  from  one  parametric  polynomial  to  the 
next.  For  a  uniform,  non-periodic  B-spline,  the  knots  are 
evenly  spaced  and  defined  as   follows: 

ti  =  0  if  i  <  k 

ti     -  i  -  k  +  1        if  k  <   i  <  n  (2.12) 

ti=n-k+2  ifi>n 

with  0   <   i  <  n  +  k 
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The  range  of  the  parameter  uis0<u<n-k+2.  in 
this  evaluation,  the  convention  0/0  =  1  is  assumed. 

The  derivatives  of  B-spline  curves  are  given  as  follows  : 
First  Derivative 


n 
P'(U)  =  S  Pi  Hi  k'(u)  (2.13) 

i=0 


where 

Ni,l'(u>  =  0  (2.14) 

„    1(u)  .   Ni,lc-l(u)    *    (u    -    tj)  Nj^^-tu) 
ti+k-1  -  ti 

(tj+k  -  U)  Nj+l,k-l'(U)  -  Nj+^fc-iJu) 

ti+k  ~   ti+1 


(2.15) 


Second  Derivative 


n 
P"(U)   =   S   Pi  Ni ,k"(u)  (2.16) 

i=0 
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where 

Ni,l"(u)    =   0  (J.17) 

Ni,2"(u)    =   0  (2.18) 

2    *    *i,k-l'(u)     H-    (u    -    tl)    Mj^-x'-Cu) 

wi,k      (u)    - 

fci+k-l  "  ti 

(tj+k    -    U)    Nj-n^-x'^u)     -    2    *    Nj+x^-x'tu) 
fci+k  "  fci+l 

(2.19) 

Third  Derivative 


n 
P"'(u)    =   2    Pi    Ni/k'"(u)  (2.20) 

1=0 


where 


Ni,l'"(u>    =    0  (2.21) 

Ni,2'"(u)    =    0  (2.22) 

Ni,3'"(U)    =    0  (2.23) 


■lfk...(U)  3    *  Ni,k-l"(u)    +    (u   -  ti)    Hi^-i'-'Cu) 


fci+k-l   "  ti 
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(tj+k    -    U)    Ni+1(]{_1'"(u)     -    3    *    Ni+1>k.1"(u) 
fci+k  "  ti+l 

(2.24) 

In  general,    the  derivative  of  degree  d  of  a  B-spline   is 

n 
Pd(u)    =  s  PjNi  kd(u)  (2# 

i=0 


where 


Ni,jd(u)    =0        ,      j    =   l,    2 d 


(2.26) 


Ni   kd(u)    .  d   *   Nj^-xCd-Dtu)    +    (u   -   ti)    Ni^.x^tu) 


fci+k-l   _  ti 

(tj+k  ~  u)    Ni+ltk_1d(u)    -  d  *  Ni+1>k_x(d-l) (u) 

fci+k    _    ti+l 

(2.27) 


2.6   B-spline  Cams 

As  discussed  in  Section  2.4,  conventional  design  meth- 
ods use  either  standard  curves  or  polynomials  to  represent 
the  cam  motion.  These  curves  are  usually  defined  in  non- 
parametric  form. 
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The  design  approach  taken  in  this  work  is  to  use  B- 
splines  to  represent  the  lift  diagram  and  to  use  this  repre- 
sentation in  the  design  of  the  cam.  Since  B-splines  are 
parametric  curves,  normalized  lift  as  well  as  the  cam  rota- 
tion now  become  functions  of  the  independent  B-spline  param- 
eter "u": 


n 


9  =  2  Pix  Ni  k(u)  (2.28) 

i=0 


y  -  S  Piy  Nj_  k(u)  (2.29) 

1=0  ' 


where 

pi,y  are  the  y-coordinates  of  the  control  points 
pi,x  are  the  9-coordinates  of  the  control  points 

Consequently,  the  true  lift  of  the  follower  is  also  a  func- 
tion of  u: 


n 
Y  =  L  *   2  Piy  Ni/k(u)  (2.30) 


In  the  design  of  conventional  cams,  the  lift  is  generally 
treated  as  an  explicit  non-parametric  function  of  cam  rota- 
tion; it  is  clear  from  the  above  equations  that  this  ap- 
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proach  has  to  be  modified  in  the  case  of  B-spline  curves. 
Here,  we  must  treat  both  lift  and  cam  rotation  as  being 
functions   of  the  parameter  u   in  all   derivations. 

In  addition  to  the  cam  angle  of  rotation  and  lift, 
their  derivatives  must  be  also  taken  into  consideration. 
The  derivatives  of  Equations  (2.28)  and  (2.30)  can  be  ex- 
pressed as   follows: 


n 
de/du   =   2   Pix   Ni/k'(u)  (2.31) 

i=0 


n 
dY/du   =    L    *    2    Piy   Ni    k'(u)  (2.32) 

1=0 


n 
d29/du2    =      S    pix   Ni/k"(u)  (2.33) 

i=0 

n 
d2y/du2    =    L    *    2   Piy   Nt    k"(u)  (2.34) 

i=0 

n 
d3e/dU3    -  S    pix   Ni(k'"(u)  (2.35) 


d3Y/du3    =    L    *    S    Piy    Ni    k'"(u)  (2.36) 

1=0 
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Reduced  velocity,  reduced  acceleration  and  reduced  jerk  are 
given  by  the  following  equations: 


dY/de  =  (dY/du)/(d9/du)  (2.37) 

(de/du)*(d2Y/du2)  -  (dY/du)*(d29/du2) 


d2Y/d62 


(de/du) 3 

(2.38) 


3    3    (de/du)* (de/du  *  d3Y/du3  -  dY/du  *  d3e/du3) 

(de/du) 5 
3*d2e/du2  (  de/du  *  d2Y/du2  -  dY/du  *  d28/du2) 


(  de/du  )5 

(2.39) 

2.7   Boundary  Conditions  for  B-spline  Cams 

In  this  work,  it  is  assumed  that  the  cam  consists  of  a 
single  rise  and  a  single  return,  with  the  boundary  condi- 
tions specified  below: 
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Rise  Zone 

9//S  =  0 

Y=O,Y'=0,Y',=O,Y,,'=O 

e//3  =  o 

Y=l,   Y1=O,Y,'=0,Y','=0 
Return  Zone 

9//3  =  0 

Y=0,Y'=0,Y,,=0,Y',,=0 

e/p   =  o 

Y=l,   Yl=0,Y,I=O,Y,l,=0 

The  remaining  of  this  section  is  devoted  to  a  discussion  on 
the  satisfaction  of  these  conditions.  It  is  to  be  noted 
that  the  implementation  of  B-splines  is  also  applicable  for 
cams  with  different  boundary  condition  specifications. 

The  approach  adopted  in  this  thesis  is  to  use  one  B- 
spline  curve  to  represent  the  rise  curve  and  one  B-spline 
curve  to  represent  the  return  curve.  The  number  of  control 
points  used  for  each  of  these  B-splines  is  11.  Some  of  the 
coordinates  of  these  points  may  be  determined  by  the  speci- 
fied boundary  conditions.  The  relationship  between  the 
boundary  conditions  and  the  control  points  is  discussed  in 
the  next  section. 
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2.7.1  Position  boundary  conditions 

As  shown  in  Figure  1.1,  a  B-spline  curve  passes  through 
both  end  control  points.  In  a  rise  zone,  position  boundary 
conditions  can  be  satisfied  by  setting  the  coordinates  of 
both  end  points  equal  to  the  coordinates  of  the  ends  of  the 
rise  cycle.  For  example,  if  the  left  end  of  the  rise  is  at 
(0,0)  and  the  right  end  is  at  (jt,1),  the  first  control 
point  must  be  set  at  (0,0),  and  the  last  control  point  at 
(tt,1).  A  similar  approach  is  applicable  for  the  return 
cycle.  In  this  way,  the  position  boundary  conditions  at  the 
end  points  of  each  cycle  are  satisfied. 


2.7.2  Derivative  boundary  conditions 

For  any  non-periodic  B-spline,  the  curve  is  tangent  to 
the  first  and  last  spans  of  the  control  polygon  at  the  first 
and  last  control  points.  This  characteristic  can  be  used  to 
satisfy  the  first  derivative  boundary  conditions.  In  order 
to  ensure  zero  velocity  boundary  conditions,  the  tangent  of 
the  curve  at  the  appropriate  boundary  must  be  set  to  zero. 
With  reference  to  Figures  2.7(a)  and  2.7(b),  this  can  be 
achieved  by  setting  the  y-coordinates  of  control  points  Pa 
and  P9  equal  to  the  y-coordinate  of  the  first  and  the  last 
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Figure  2.7(a)  :  a  Rise  Cycle  B-spline 


P7P8P9P10 
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Figure  2.7(b)  :  A  Return  Cycle  B-spline 
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points  respectively.  A  similar  approach  can  be  applied 
to  acceleration  and  jerk  boundary  conditions.  To  obtain 
zero  acceleration  and  jerk  at  both  boundaries,  the  y-coordi- 
nates  of  points  P2  and  P3  must  be  set  equal  to  the  y-coordi- 
nate  of  point  P0,  while  the  y-coordinates  of  points  P7  and 
P8  must  be  set  equal  to  y-coordinate  of  point  Pio-  As  a 
result  of  imposing  the  boundary  conditions,  the  x-coordi- 
nates  of  points  P0  and  P10,  and  the  y-coordinates  of  points 
pl<  p2'  p3>  p7'  p8/  and  p9  are  fixed.  The  x-coordinates  of 
points  plf  P2,  P3,  P4,  P5,  p6/  p7,  p8  and  P9 ,  and  the  y- 
coordinates  of  points  P4 ,  P5  and  P6  are  free  variables. 
They  are  used  as  design  variables  to  satisfy  constraints  and 
to  reduce  cost  in  an  optimization  process. 

2.8  Advantages  of  B- spline  Curves 

B-spline  curves  possess  several  properties  which  make 
them  very  attractive  for  cam  design  applications. 

a)  Degree  of  continuity  :  The  degree  of  continuity  of  a 
B-spline  curve  is  controlled  only  by  the  parameter  k. 
With  a  specified  value  of  k,  the  continuity  of  the 
curve  is  (k-2) .  Therefore,  the  desired  degree  of 
continuity  can  be  obtained  by  setting  k  to  an  appro- 
priate value.  For  example,  if  continuity  up  to  jerk 
is  required,  the  value  of  k  must  be  at  least  5.  In 
this  work,  the  value  of  k  is  set  to  5. 
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b)  Flexible  number  of  control  points:  As  noted  earlier, 
the  number  of  control  points  used  in  this  work  is  11. 
The  boundary  conditions  impose  restrictions  on  the 
coordinates  of  some  of  the  control  points,  while  the 
remaining  coordinates  can  be  treated  as  design  varia- 
bles. B-spline  curves  can  be  extended  to  have  more 
control  points  without  affecting  local  properties.  As 
a  result,  more  points  can  be  added  to  provide  more 
design  variables  to  satisfy  more  rigorous  performance 
requirements . 

c)  Good  local  control:  B-splines  provide  good  local 
control.  Local  control  is  the  ability  to  change  a 
certain  section  of  a  curve  without  affecting  the  shape 
of  the  curve  far  away  from  the  region  of  change.  As 
shown  in  Figure  2.8,  the  location  of  control  point  P5 
is  shifted  to  a  new  position  at  P5*.  The  effect  of 
shifting  the  point  changes  only  a  small  neighboring 
section  of  the  curve.  This  indicates  that  B-spline 
curves  provide  good  local  control.  Such  control  plays 
an  important  role  in  the  design  of  cams.  Very  often, 
only  a  few  sections  of  the  profile  of  the  initial 
design  violate  the  design  requirements.  It  is  there- 
fore desirable  to  retain  the  sections  with  no  viola- 
tion, and  correct  only  the  problem  regions.   This  is 
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Figure  2.8:   Illustration  of  Local  Control 
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not  possible  if  either  standard  curves  or  polynomials 
are  used,  since  they  do  not  provide  local  control.  By 
using  B-splines,  the  violation  may  be  reduced  or 
eliminated  by  shifting  the  control  points  near  the 
problem  regions  appropriately,  while  retaining  the 
shape  of  the  curve  at  other  regions  where  no  violation 
occurs . 
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CHAPTER  III 
DESIGN  CONSIDERATIONS 

For  a  cam  design  to  be  acceptable,  it  must  satisfy  a 
specified  set  of  performance  requirements.  Therefore,  it  is 
necessary  to  establish  a  set  of  design  relationships  to 
account  for  these  requirements.  These  performance  require- 
ments are  enforced  in  the  optimization  process  as  con- 
straints. The  design  requirements  considered  in  this  thesis 
are  pressure  angle,  separation,  oscillation,  stress  and 
undercutting.  The  details  of  these  considerations  are 
discussed  below.  Other  design  considerations  are  treated  as 
problem  dependent,  with  the  designer  being  allowed  to  add  as 
many  other  performance  constraints  as  required. 

3.1   Pressure  Angle 

The  pressure  angle  is  the  angle  between  the  normal  to 
the  pitch  curve  and  the  axis  of  the  follower  stem.  The 
value  of  the  pressure  angle  should  be  kept  low  to  reduce 
friction  and  chattering  of  the  follower,  as  well  as  to 
obtain  better  force  transmission  [1].  From  Figure  2.3,  we 
observed  that  only  the  force  component  along  the  axis  of  the 
follower  is  used  to  overcome  the  output  load.  The  tangen- 
tial component  of  the  force  causes  friction  between  the 
follower  stem  and  its  guideway.   Therefore,  a  smaller  pres- 
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sure  angle  provides  better  force  transmission,  and  increases 
the  smoothness  of  the  follower  motion. 

Designers  often  limit  the  pressure  angle  to  a  maximum 
of  30°.  However,  if  the  follower  bearings  are  good  and  the 
cam  follower  is  rigid,  the  maximum  pressure  angle  may  be 
increased  to  a  larger  value. 

In  Figure  2.3,  using  the  Aronhold-Kennedy  theorem  of 
three  centers,  the  instantaneous  center  of  velocity  between 
cam  2  and  follower  4  is  at  P24.  since  the  follower  is 
translating,  all  points  of  the  follower  have  velocity  equal 
to  that  of  P24,  which  is  also  equal  to  the  velocity  of  the 
coincident  point  on  link  2  [3]. 

VP  =  Y  =u  *  Rp  (3-1) 

where 

w     is  the  angular  velocity  of  cam 

Rp     is     the     distance     between     cam     center     of 
rotation   and  point   P24 

Dividing   equation    (3.1)    by w  ,    the   equation   is   reduced   to 

Y/w  =   Y"    =   Rp   =   e   +    (a   +   Y)    tan(<£)  (3.2) 

where  a  is  the  vertical  distance  from  the  cam  axis  to  the 
prime  circle  and  is  given  by 
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a  =  sqrt(  RQ2  -  e2)  (3-3) 

where 

R0  is  base  radius  of  the  cam 

e  is  the  follower  offset. 

The  final  expression  for  pressure  angle  is  given  by  : 

Y"  -  e 

*  "   tan  _1  . (3.4) 

sqrt(  RQ2  -  e2  )+  V 


3.2   Separation 

The  cam  and  the  roller  follower  must  maintain 
continuous  contact  at  all  times.  Failure  to  maintain 
continuous  contact  between  the  cam  and  the  follower  will 
result  in  separation.  This  is  generally  undesirable  since 
the  motion  generated  by  the  roller  follower  may  be  unaccept- 
able, in  addition,  separation  may  cause  other  problems  such 
as  impact  between  the  follower  and  cam  surface,  wear,  noise, 
etc.  To  avoid  separation,  we  examine  the  dynamics  of  the 
cam  system.  As  shown  in  Figure  2.3,  the  driving  force  of 
the  cam  is  given  by 

Fd  "  (  PL  +  Mf  *  Y  )/cos(0)  (3.5) 
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where 

PL  is  the  preload  specified  by  the  user  (including 
friction,  spring  force,  etc.) 

Mf  is  the  mass  of  the  follower  mechanism 
V  is  the  true  acceleration  of  the  follower 

The  force  F<j  must  always  maintain  a  positive  value 
for  continuous  contact  between  the  cam  and  the  follower. 
However,  if  the  acceleration  required  to  track  the  cam  is 
too  high  on  the  return  cycle,  a  negative  value  of  F,j  is  re- 
quired to  maintain  contact.  Since  this  is  physically  impos- 
sible, it  will  result  in  separation  between  the  cam  surface 
and  the  follower. 

To  prevent  separation,  one  option  is  to  increase  the 
preload  [2].  However,  an  increase  in  the  preload  will 
increase  the  contact  stress  between  the  cam  and  the  follow- 
er. A  second  alternative  is  to  keep  the  acceleration  during 
the  return  cycle  to  a  minimum  by  controlling  the  shape  of 
the  profile.  In  this  thesis,  separation  between  the  follow- 
er and  the  cam  is  avoided  by  imposing  a  constraint  in  the 
optimization  process,  which  restricts  the  driving  force  to 
be  positive  throughout  the  motion. 
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3.3   Elimination  of  Oscillation  in  the  Follower  Velocity 

Oscillation  is  a  condition  in  which  the  follower  at- 
tempts to  descend  during  the  rise  cycle,  or  ascend  during 
the  return  cycle.  Due  to  oscillation,  it  is  possible  that 
in  the  interior  of  a  rise  or  a  return,  the  follower  may 
attempt  to  move  to  a  location  beyond  the  specified  lower  or 
upper  boundary  limits.  This  is  unacceptable  and  should  be 
avoided,  since  the  cam  system  may  be  confined  within  a 
limited  space  and  movements  beyond  the  limits  may  cause 
severe  damage  to  the  system.  Therefore,  it  is  necessary  to 
design  a  cam  with  monotonic  ascent  during  the  rise  and 
monotonic  descent  during  the  return  to  avoid  oscillation. 

In  this  work,  monotonic  ascent  during  the  rise  and 
monotonic  descent  during  the  return  are  obtained  by  re- 
stricting the  velocity  of  the  cam  follower  to  be  positive 
throughout  the  rise  cycle  and  negative  throughout  the  return 
cycle. 

3.4   Stress 

Stress  at  the  area  of  contact  between  the  cam  and  the 
follower  is  another  important  factor.  Excessive  stress  can 
cause  surface  wear  and  failure  of  the  cam  system. 

According  to  Hertz  [1],  the  stress  distribution  of  two 
cylinders  in  contact,  assuming  perfect  alignment,  is  rectan- 
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gular  in  shape.  The  maximum  compressive  stress  of  two 
perfectly  aligned  cylinders  is  at  the  center  of  the  contact 
area  and  decreases  to  zero  at  the  ends. 

The  Hertzian  equation  to  calculate  the  maximum  compres- 
sive stress  for  two  cylinders  in  contact  and  in  alignment 
is  given  as  follows  [1]: 


.  ,  _   Fd   (  rc     +  rf  ) 

SC   -         ■ (3.6) 

7r*tc*rf  *rc*  (  ( 1-MC2 )  /Ec+  ( 1-Mf 2 )  /Ef ) 


where 


Sc  =  Maximum  compressive  stress  at  any  point,  lb/ in.2 

Fd  =  normal  load  on  cam  profile,  lb. 

rc  and   rf   =  radii   of  curvature   of   cam   and 

follower  respectively,  in. 
tc  =  thickness  of  contacting  cam  and  follower 
Mc  and  Mf  =  Poisson's  ratio   for  cam  and   follower 

respectively. 
Ec  and   Ef  =  moduli   of  elasticity  of   cam   and 
follower,  respectively. 
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3 . 5   Undercutting 

The  radius  of  curvature  of  the  cam  profile  is  another 
important  consideration  in  the  design  of  the  cam.  An  ade- 
quate cam  curvature  is  required  to  produce  correct  follower 
movement.  if  the  curvature  of  the  cam  is  too  sharp,  the 
follower  may  not  be  able  to  trace  the  profile  correctly,  and 
thus  the  motion  generated  may  be  unacceptable.  In  addition, 
a  sharp  curvature  may  cause  excessive  stress  on  the  cam  sur- 
face, which  may  cause  wear  and  eventual  failure  of  the 
system.   This  condition  is  called  undercutting. 

With  a  roller  follower,  undercutting  occurs  in  a  convex 
curve  when  the  radius  of  curvature  of  the  pitch  curve  is 
less  than  the  radius  of  the  roller.  With  a  concave  curve,  it 
occurs  when  the  radius  of  curvature  of  the  pitch  curve  is 
less  than  zero  or  the  radius  of  curvature  of  the  cam  profile 
is  less  than  the  radius  of  the  roller.  Figure  3.1  illus- 
trates the  problems  of  undercutting.  For  a  safe  design,  the 
curvature  of  the  profile  must  be  greater  than  one  and  a  half 
times  the  radius  of  the  follower  [9].  The  equation  for  the 
radius  of  curvature  is  given  as  follow  : 


(  c2  +  (yi  -  e)2  )  3/2 


(3.7) 


c2  +  (Y1  -  e)  (2*V  -  e)  -  c*Y" 
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No  Undercutting 


Figure  3.1  :   Illustration  of  Undercutting 

(  Extracted  from  Reference  [14]  ) 
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where 

c  is  the  distance  between  cam  center  and  cam  profile, 
given  by 

c  =  a  +  X  (3.8) 

Undercutting  can  be  avoided  by  either  or  both  of  the  follow- 
ing ways  : 

1)  Use  a  small  roller.   However,  a  small   roller  size 
may  cause  excessive  stress  in  the  roller. 

2)  Use  a  larger  cam.   This  approach  generally  has 

a   practical   limitation,   since  a   larger   cam 
increases  space  requirements,  material  cost,  etc. 

In  this  thesis,  the  radius  of  the  follower  is  kept  constant; 
therefore  undercutting  is  to  be  avoided  by  properly  sizing 
the  cam. 

The  considerations  that  have  been  described  are  en- 
forced as  constraints  in  the  optimization  process  as  de- 
scribed in  Chapter  V. 
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CHAPTER  IV 
MANUFACTURING  ASPECTS 


In  this  research,  our  purpose  is  to  simultaneously 
consider  the  design  and  manufacturing  aspects  in  order  to 
design  cams  with  the  lowest  manufactured  cost  or  lowest 
cutting  time.  For  designing  cams  with  low  machining  time,  a 
mathematical  model  for  estimating  the  cutting  time  must 
first  be  derived.  The  relevant  equations  for  estimating  the 
cutting  time  are  discussed  later  in  this  chapter.  For 
minimizing  the  manufactured  cost,  we  consider  the  cost  of 
the  metal  required  to  make  the  cam  to  be  the  most  signifi- 
cant contributor  to  the  total  cost.  Since  cams  with  small- 
er size  require  less  material,  it  follows  that  a  minimum 
size  cam  minimizes  the  total  cost. 


4.1  cutting  Time 

There  are  a  variety  of  ways  to  manufacture  plate  cams. 
One  possibility  is  the  use  of  numerically  controlled  milling 
machines,  which  offer  significant  benefits  in  terms  of 
speed,  reliability  and  accuracy.  The  use  of  NC  machines 
requires  the  process  planner  to  generate  data  points  for  the 
tool  path.  Today,  with  the  aid  of  digital  computers,  the 
data  points  for  the  tool  path  can  be  generated  conveniently, 
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and  then  transferred  to  a  tape  and  fed  into  the  NC  machine. 
In  the  present  work,  we  consider  manufacture  of  the  cam 
using  a  Pratt  &  Whitney  Tapemate  Series  C  NC  milling  ma- 
chine. This  is  a  bed-type  vertical  milling,  drilling  and 
tapping  machine  with  a  7-tool  turret.  The  machining  method 
selected  to  manufacture  the  cam  is  the  end  milling  slotting 
process.  Once  the  geometry  of  the  cam  is  decided  by  the 
optimization  process,  a  blank  is  then  cut  to  the  desired 
geometry  and  shape. 

During  the  manufacturing  process,  it  is  necessary  to 
select  a  set  of  machining  parameters,  such  as  tool  size, 
spindle  speed,  etc.  These  parameters  are  generally  decided 
based  on  the  engineer's  experience.  For  the  purpose  of  this 
work,  the  selection  of  the  necessary  machining  parameters 
are  made  by  referring  to  standard  tables  [10].  As  indicated 
in  Table  1,  a  tool  diameter  of  3/4  inches  provides  a  good 
feed  rate  for  a  wide  variety  of  materials.  Furthermore,  a 
tool  size  of  3/4  inches  is  also  the  most  appropriate  consid- 
ering both  tool  life  and  productivity. 

The  other  relevant  manufacturing  parameters  that  were 
selected  are  listed  below: 

Tool  material  High  Speed  Steel 

Tool  diameter 3/4  inches 

Number  of  teeth  4 
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Table  1  :   Machining  Data  for  End  Milling  Proces 


MATERIAL 
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COT 
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(  Extracted  from  Reference  [10]  ) 
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Depth  of  cut 0.25  inches 

Feed  per  tooth  0.002  inches 

Cutting  speed no  feet  per  minute 

Once  the  cutting  speed  and  tool  diameter  are  known,  the 
spindle  speed  in  RPM  can  be  calculated  as  follows  [11]: 

Vs  =  (12  *  Vc)  /  (*  *  dt)  (4.1) 

where 

Vc  is  the  nominal  cutting  speed  in  fpm 
dt  is  the  tool  diameter  in  inches 

The  Pratt  &  Whitney  machine  used  in  this  work  has  only 
six  spindle  speeds.  Therefore,  the  closest  lower  speed  is 
selected.  Using  this  spindle  speed,  the  actual  cutting 
speed  is  given  by 

Vca  =  12  /  (60  *  w    *   dt)  *  Vs'  (4.2) 

where 

Vca  is  the  actual  cutting  speed  in  fpm 
Vs '  is  the  actual  spindle  speed  in  RPM 
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Now,  the  feed  in  inches  per  second  can  be  computed  as 
follows: 

Feed  =  fpt  *  Nt  *  Vs'/60  (4.3) 

where 

fpt  is  the  feed  per  tooth  in  inches 
Nt  is  the  number  of  teeth  in  the  tool 

4.1.1  Generation  of  data  points 

The  machining  of  the  cam  requires  the  computation  of 
data  points  for  the  tool  path.  Conventionally,  a  series  of 
fixed  steps  is  taken  along  the  x-axis  and  y-axis  of  the 
milling  machine  to  generate  the  cam  contour  [2].  The  size 
of  the  steps  is  determined  at  the  most  critical  section  of 
the  profile,  where  the  step  size  is  required  to  be  suffi- 
ciently small  to  meet  the  desired  accuracy.  However,  it  is 
possible  to  use  variable  step  sizes  to  generate  the  contour, 
because  at  some  sections  a  larger  step  may  be  allowed  for 
the  same  accuracy  [4].  The  advantage  of  this  approach  is 
that  it  minimizes  the  number  of  turns  of  the  tool  path,  thus 
reducing  the  cutting  time  significantly. 

As  shown  in  Figure  4.2,  the  data  points  are  restricted 
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s=.001 


Desired  Surface 

Possible  Stepped  Surface 

•Stepping  Increment 


esi red  Curve 


Figure  4.1  :  Fixed  Step  Size  Path  Generation 
(  Extracted  from  Reference  [2]  ) 
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Figure  4.2  :  Variable  Step  size  Path  Generation 
(  Extracted  from  Reference  [4]  ) 
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to  be  on  the  outside  of  the  cam  profile.  Therefore,  the 
cutter  never  cuts  into  the  profile,  but  remains  outside 
within  a  specified  tolerance  "e".  This  strategy  will  allow 
subsequent  improvement  of  the  surface  finish  and  accuracy  of 
the  cam  by  passing  it  through  a  finishing  process,  such  as 
profile  grinding. 

The  cam  motion  consists  of  a  rise,  a  top  dwell,  a 
return  and  a  bottom  dwell.  The  approach  adopted  is  to 
generate  data  points  for  each  cycle  individually.  There- 
fore, the  starting  point  of  the  tool  path  corresponding  to  a 
cycle  is  at  the  start  of  the  cycle,  and  the  ending  point  of 
the  tool  path  is  at  the  other  end  of  the  cycle. 

To  generate  data  points,  it  is  necessary  to  divide  the 
profile  into  different  sections,  based  on  the  sign  of  the 
slope  dY/dX.  This  is  done  by  computing  the  critical  points 
of  the  profile  with  respect  to  the  coordinate  axes,  i.e., 
points  at  which  the  slope  is  either  zero  or  infinity.  Each 
section  is  confined  between  two  consecutive  critical  points 
as  shown  in  Figure  4.2.  Note  that  the  sign  of  the  slope 
remains  the  same  within  a  section. 

The  purpose  of  dividing  the  cam  profile  into  different 
sections  is  because  the  sequence  of  data  point  generation 
for  a  section  with  positive  slope  is  different  from  a  sec- 
tion with  negative  slope: 
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a)  Figure  4.3(a)  shows  a  section  with  positive  slope. 
Starting  at  point  (X0,Y0),  the  tool  moves  along  the  x- 
axis  to  point  (XlfYx),  so  that  the  maximum  deviation 
from  the  profile  is  equal  to  the  specified  tolerance 
"e".  The  tool  then  moves  along  the  y-axis  to  (X2,Y2). 
Notice  that  this  strategy  maximizes  the  step  sizes  along 
each  direction,  and  thus  reduces  the  number  of  turns  in 
the  tool  path.  This  procedure  is  repeated  until  the 
entire  section  is  completed. 

b)  For  a  section  with  negative  slope  as  shown  in  Figure 
4.3(b),  starting  from  (X0,Y0),  the  tool  first  moves 
along  the  y-axis  to  (X1,Y1) ,  and  then  along  the  x-axis 
to  (X2,Y2). 

In  computing  the  data  points,  since  the  rise  and  the 
return  are  represented  by  parametric  curves,  the  process  of 
data  point  generation  for  these  cycles  is  different  from  the 
process  that  is  used  for  the  dwells. 

1)  Rise  and  Return 

The  process  of  data  point  generation  is  executed  as 
follows: 
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Figure  4.3(a)  :  Generation  of  Data  Points  for  a  section 
with  Positive  Slope 
(  Extracted  from  Reference  [4]  ) 
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Figure  4.3(b)  :  Generation  of  Data  Points  for  a  Section 
with  Negative  slope 
(  Extracted  from  Reference  [4]  ) 
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a)   If  dY/dX  >  0  : 

Figure  4.3(a)  shows  a  section  of  profile  with  positive 
slope.  The  geometric  relationship  between  points  (X0,Y0) 
and  (Xi,Y]J  is  given  by 

f(u)   =  (  r(u)  +  e    )  *  sin(  9(u)  )  -  Y0  =   0       (4.4) 

where 

r(u)   is   the  instantaneous  radial  distance  to  the   cam 

profile 
8(u)  is  the  angle  corresponding  to  the  data  point 

The  equations  for  r(u)  and  9(u)  are  given  as  follows: 


n 
r(u)  =  R0  +  L*  S  PiyNi  k(u)  (4.5) 

i=o 


n 
8(u)  =  S  Pix  Ni  k(u)  (4.6) 

1=0  ' 


Similarly,  a  mathematical  expression  can  be  derived  from  the 
geometric  relationship  between  points  (X^YJ  and  (X2,Y2): 

g(u)  =  r(u)  *  cos(  e(u)  )  -  X1  =   o  (4.7) 
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Referring  to  Figure  4.3(a),  suppose  that  X0  and  Y0  are 
known  coordinates,  and  the  angle  at  this  point  is  also 
known.  To  generate  the  next  data  point  (Xl,Y1) ,  a  Newton- 
Raphson  iteration  is  implemented  to  solve  Equation  (4.4). 

uj+1  -  Uj  -  f(Uj)/f  (Uj)  (4.8) 

where 

j  is  the  iteration  number 

and  f'(Uj)  is  given  by 

f'(uj)  =  (  r(Uj)  +  6  )  *  cos(  e(Uj)  )  *  e'(uj) 

+  r'(Uj)  *  sin(  8(uj)  )  (4.9) 


where 


n 
6'(Uj)  =  Z   Pix  Nirk'(Uj)  (4.10) 

1=0 


When  the  parametric  value  u  is  known,  the  coordinates  of 
(Xi,Yi)  are  given  by 

Xj.  =  r(u)  *  cos(  9(u)  )  (4.11) 

Yl  "  Y0  (4.12) 
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The  next  data  point  to  be  generated  is  (X2,Y2).  The 
value  of  the  parameter  u  corresponding  to  this  point  can  be 
calculated  by  the  following  Newton-Raphson  iteration  for 
solving  Equation  (4.7) 

uj+1  =  uj  -  g(Uj)/g'(uj)  (4.13) 

where 

j  is  the  iteration  number 

and  g'(uj)  is  given  by 

g' (Uj)  =  r' (uj)  *  cos(  6(uj)  ) 

-  r(uj)  *  sin(  9(uj)  )  *  9'(Uj) 

(4.14) 

Once  the  value  of  u  is  known,  the  coordinates  of  point 
(X2<Y2)  are  given  by 

X2  =  xl  (4.15) 

Y2  =  r(u)  *  sin(  9(u)  )  (4.16) 

b)   If  dY/dX  <  o  : 

Figure    4.3(b)     shows    a    section    of    cam    profile    with 
positive   slope.      The   geometric   relationship   between   points 
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(X0,Y0)  and  (X^Yi)  can  be  expressed  as 

f(u)  =  (  r(u)  +  6  )  *  cos(  6(u)  )  -  X0  =  0  (4.17) 
Points  (Xi.Yx)    and  (X2,Y2)  are  related  as  follow: 

g(u)  =  r(u)  *  sin(  6(u)  )  -  Yx  =  0  (4.18) 

As  was  the  case  for  a  section  with  positive  slope,  data 
points  for  a  section  with  negative  slope  are  generated  using 
a  Newton-Raphson  iteration  to  solve  Equations  (4.17)  and 
(4.18). 

The  algorithm  for  tool  path  generation  was  used  to 
generate  NC  code  automatically.  The  code  was  then  trans- 
ferred to  a  tape  and  fed  into  the  Pratt  &  Whitney  Tapemate 
Series  C  NC  machine  to  produce  the  cam. 

2)  Dwells 

The  function  needed  to  represent  dwell  cycles  is  simple 
and  in  non-parametric  form.  In  this  thesis,  the  procedure 
for  data  point  generation  for  dwell  cycles  is  taken 
directly  from   [4] . 

4.1.2     Estimation  of  cutting  time 

In  machining  the  cam  profile,  the  number  of  turns  in 
the  tool  path  is  a  significant  factor  in  estimating  the 
cutting   time.      A   series   of   tests   was   conducted   on   the   Pratt 
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&  Whitney  NC  machine  to  estimate  the  time  required  to  make 
a  turn  [4].  The  NC  machine  was  programmed  to  travel  along 
different  paths  of  equal  length  but  with  different  numbers 
of  turns.  These  tests,  which  are  fully  described  in  [4], 
are  summarized  below. 

The  first  path  selected  was  along  the  sides  of  a  square 
of  side  4  inches  as  shown  in  Figure  4.4.  The  time  taken  by 
the  tool  to  traverse  the  sides  of  the  square  for  a  total  of 
20  cycles  was  measured.  Note  that  the  tool  is  required  to 
make  4  turns  per  cycle  on  this  path  and  the  length  of  the 
path  is  16  inches  per  cycle. 

In  the  next  test,  the  tool  was  driven  along  a  staircase 
path  from  corner  1  to  corner  3,  and  back  to  corner  1  of  the 
square,  as  shown  in  Figure  4.4.  The  tool  was  thus  pro- 
grammed to  make  20  turns  per  cycle  along  this  path,  which 
also  has  a  length  of  16  inches  per  cycle.  The  total  time 
required  to  travel  2  0  cycles  on  this  path  was  also  measured. 

In  subsequent  tests,  the  tool  was  programmed  to  travel 
along  the  similar  staircase  paths,  but  with  40,  80  and  100 
turns  per  cycle.  The  distance  per  cycle  for  all  these  paths 
was  kept  equal  (16  inches  per  cycle),  even  though  the  number 
of  turns  for  each  path  varied.  The  data  recorded  is  repro- 
duced in  Table  2.  Based  on  these  results,  the  estimated 
time  to  make  a  turn  is  taken  to  be  0.4  seconds. 
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Figure  4.4  :  Estimation  of  Time  Lapse  Per  Turn  in  Tool  Path 
(  Extracted  from  Reference  [4]  ) 
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Table  2  :  Estimate  of  Time  Lapse  Per  Turn 


Case   1 


No.  of  cycles       No.  cf  tums/cycl    Total  time  (Min) 


20 
20 


14.133 
14.155 


Case 


Number  of   Number  of     Total     Total  time    Time 
cycles     turns/cycl    time (min)   lapse (min)  lapse/turn 


20 

20 

20 

40 

20 

80 

20 

100 

16-35  2.216  .33seconds 

19.35  5.216  .39seconds 

23.947  9.81  .37seconds 

39.633  25.58  .38seconds 


(  Extracted  form  Reference  [4]  ) 
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Once  the  path  generation  has  been  completed,  the  dis- 
tance the  tool  needs  to  travel  to  cut  the  cam  profile  can  be 
computed  as  follows: 


n 

Sd  =  tl0l    (Yi  "  Yi"l)2  +  (Xi  "  Xi"l>2  11/2       <4-19> 

where 

Xi  is  the  x-coordinate  of  the  data  point  "i" 
Yi  is  the  y-coordinate  of  the  data  point  "i" 
n  is  the  total  number  of  steps 

Now,  the  total  time  required  to  cut  the  cam  can  be  estimated 
as 


Tcut  =  Sd/Feed  +  0.4  *  Nt  (4.20) 

where 

Nt  is  the  total  number  of  turns 

Equation  (4.20)  is  used  as  a  cost  function  in  the  optimiza- 
tion process. 
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4.2   Cross-sectional  Area  of  Cam 

The  material  required  to  manufacture  a  cam  is  directly 
related  to  the  cam  volume.  Assuming  constant  depth, 
the  minimization  of  cam  material  is  equivalent  to  the  mini- 
mization of  cam  area. 

In  a  previous  chapter,  we  mentioned  that  the  cam  pro- 
files of  interest  consist  of  a  rise,  a  return,  a  top  dwell 
and  a  bottom  dwell.  Therefore,  the  total  cam  area  can  be 
obtained  by  adding  up  the  area  of  these  zones. 


Al  +  A2  +  A3  +  A4 


(4.21) 


where 


A  is  the  total  cam  area 

Ai  is  the  area  for  the  rise 

A2  is  the  area  for  the  return 

A3  is  the  area  for  the  top  dwell 

A4  is  the  area  for  the  bottom  dwell 

Equations  for  A1(  A2,  A3  and  A4  are  given  as  follows: 


0.5  *  r2(9)  d9  (4.22) 

67 


A2= 


73 1 
0.5  *  r2(6)  de 


(4.23) 


A3=  0.5  *  (RQ  +  L)2  *  et 
A4=  0.5  *  RQ2  *  eb 
where 


(4.24) 
(4.25) 


r(9)  is  the  instantaneous  radial  distance  to  the 

cam  profile 
R0  is  the  base  radius  of  cam 
L  is  the  total  lift 

«1,  Pi   are  the  bounds  for  rise  and  return  zones 
9-t  is  the  total  angle  for  top  dwell 
6t,  is  the  total  angle  for  bottom  dwell 

Converting  equations  (4.22)  and  (4.23)  into  the  B-Spline 
form,  we  have 


Al  = 


0.5  *  r2(u)  *  9'(u)  du 


(4.26) 
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*2  = 


"1 

0.5  *  r2(u)  *  9' (u)  du 


(4.27) 


where 

r(u)  is  the  instantaneous  radial  distance  to  the  cam 

profile,  given  by  Equation  (4.5) 
9'(u)  is  the  derivative  of  the  angle  of  rotation,  given 
by  Equation  (4.10) 

The  integrals  in  Equations  (4.26)  and  (4.27)  are  evalu- 
ated numerically  using  Simpson's  rule. 
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CHAPTER  V 
OPTIMIZATION 

In  engineering  design  applications,  the  goal  of  the 
designer  is  to  produce  the  "best"  design  of  the  desired 
system  or  component.  Numerical  optimization  is  one  of  the 
tools  that  can  be  used  to  aid  the  designer  in  this  endeavor. 

In  order  to  use  optimization  to  solve  a  given  problem, 
the  problem  must  be  stated  in  the  following  standard  format: 

Suppose  the  design  of  a  system  is  specified  by  a  vector 
b  of  r  design  variables: 

b  =  [  bx,  b2,  b3 br]T 

Find  the  value  of  b  to  minimize  a  given  cost  function,  F(b) , 
subject  to  the  constraints 

gj(b)  <  0  j  =  l,m 

hi(b)  =o  1  =  l,n 

The  optimization  problem  is  fully  defined  once  the  set 
of  design  variables,  cost  and  constraint  functions  are 
specified. 

The  optimal  design  procedure  for  solving  this  standard 
problem  is  by  iterative  improvement  of  the  design  vector,  b. 
Each  iteration  of  this  solution  process  consists  of  three 
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distinct  steps: 

1)  Analysis:  In  this  step,  the  behavior  of  the  system  is 
analyzed  at  the  current  design.  The  cost  and  constraint 
functions  are  also  evaluated  at  this  step. 

2)  Design  sensitivity  analysis:   The  derivatives  of  the 
cost  and  constraint  functions  with  respect  to  the  design 
variables  are  computed  at  this  step. 

3)  Design  update:   The  cost  and  constraint  function  values 
from  the  system  analysis  are  supplied  to  an  optimization 
algorithm  along  with  the  design  sensitivity  information. 
The  algorithm  computes  the  required  change  in  design  and 
returns  the  updated  design  for  the  next  iteration. 

Implementing  optimization  techniques  in  cam  design 
possesses  several  advantages  over  the  conventional  design 
approach.  In  optimization,  once  the  standardized  problem  is 
formulated,  they  are  several  readily  available  nonlinear 
programming  algorithms  for  solving  the  problem  automatically 
[6,  7].  On  the  other  hand,  a  conventional  design  approach 
is  tedious  and  often  requires  a  lot  of  repetitive  work.  As  a 
result,  implementing  optimization  techniques  in  cam  design 
not  only  reduces  the  amount  of  designer  time  needed,  but 
also  reduces  the  level  of  skill  required  to  produce  the 
design . 
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In  addition,  optimization  allows  integration  of  design 
and  manufacture.  This  will  ensure  that  the  final  design 
satisfies  both  design  and  manufacturing  requirements. 

5.1  Choice  of  Design  Variables 

The  first  step  in  setting  up  an  optimization  model  for 
a  design  problem  is  the  selection  of  a  set  of  design  varia- 
bles. The  designer  must  carefully  determine  which  parame- 
ters should  be  treated  as  design  variables  and  which  param- 
eters should  be  kept  constant.  If  the  number  of  design 
variables  is  insufficient,  the  cost  reduction  may  be  limit- 
ed, and  in  some  cases,  an  acceptable  design  may  be  impossi- 
ble to  obtain.  However,  if  the  number  of  design  variables 
is  large,  it  will  result  in  a  more  complex  problem.  The 
computer  time  required  will  also  increase  dramatically  and 
the  design  problem  may  be  too  expensive  to  optimize. 

In  the  formulation  used  in  this  thesis,  we  have  select- 
ed the  cam  base  radius,  the  orientation  of  the  cam  on  the 
table  of  the  NC  machine,  the  offset  of  the  roller  follower 
and  the  coordinates  of  the  control  points  of  the  rise  and 
return  B-splines  as  the  design  variables.  These  design 
variables  are  explained  in  more  detail  in  the  following 
sections. 
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1)   The  cam  base  radius 

The  cam  base  radius  determines  the  size  of  cam.  It  is 
generally  desirable  to  design  a  cam  to  be  as  small  as  possi- 
ble. A  small  cam  has  several  advantages;  it  occupies  less 
space,  it  is  cheaper  to  produce,  and  it  reduces  surface  wear 
and  torque  [1].  In  contrast,  a  large  cam  requires  more 
space  and  is  more  expensive  to  produce. 

However,  if  the  size  of  the  cam  is  too  small,  it  may 
result  in  large  pressure  angles,  high  contact  forces  and 
stresses,  and  undercutting. 

From  the  above  discussion,  it  is  obvious  that  the  base 
radius  of  the  cam  should  be  considered  as  a  design  variable. 

2)   Offset  Of  The  Roller  Follower 

From  Equation  (3.4),  it  is  observed  that  the  pressure 
angle  is  a  function  of  the  follower  offset.  During  the 
rise,  when  the  velocity  y  is  positive,  increasing  offset 
reduces  the  pressure  angle.  However,  increasing  offset 
increases  the  pressure  angle  during  the  return,  when  Y'  is 
negative.  On  the  other  hand,  decreasing  offset  reduces  the 
pressure  angle  on  the  return,  but  at  the  expense  of  an 
increased  pressure  angle  during  the  rise.  We  consider 
offset  as  a  design  variable  which  is  used  to  satisfy  pres- 
sure angle  requirements. 


73 


3)   Orientation  Of  The  Cam  Blank 

In  Chapter  IV,  we  discussed  the  procedure  for  variable 
step  size  path  generation  implemented  in  this  thesis.  The 
tool  path  is  generated  by  alternately  moving  the  tool  along 
the  x-axis  and  y-axis  of  the  NC  machine  by  varying  dis- 
tances. Since  the  steps  generated  for  the  cutter  are  along 
the  machine  axes,  the  pattern  of  steps  generated  is  directly 
dependent  on  the  orientation  of  the  cam  with  respect  to  the 
machine  axes.  This  orientation  is  measured  as  the  angle 
between  the  x-axis  of  the  NC  machine  and  the  line  connecting 
the  cam  center  of  rotation  and  the  starting  point  of  the 
rise.  For  a  fixed  profile,  varying  this  angle  may  change 
the  total  number  of  turns  and  the  total  length  of  the  cutter 
path,  since  these  are  key  factors  in  determining  the  cut- 
ting time,  the  orientation  of  the  cam  blank  will  signifi- 
cantly influence  the  cutting  time.  Therefore,  the  angle 
defining  the  orientation  of  the  blank  is  treated  as  a  design 
variable. 

Note  that  the  orientation  of  the  coordinates  axis  does 
not  affect  the  feasibility  of  the  design  with  respect  to 
performance  constraints.  Also,  this  design  variable  is  only 
applicable  for  the  minimization  of  cam  cutting  time  and  not 
for  the  minimization  of  cam  area. 
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4)   The  coordinates  of  the  control  points 

As  discussed  in  Chapter  II,  the  shape  of  the  cam  pro- 
file is  defined  by  a  series  of  control  points.  The  coordi- 
nates of  the  control  points  are  treated  as  design  variables 
in  the  optimization  for  the  satisfaction  of  constraints  and 
the  reduction  of  cost.  However,  not  all  the  coordinates  can 
be  used  as  design  variables.  Some  of  the  coordinates  may  be 
determined  uniquely  by  the  specified  boundary  conditions. 
For  example,  the  position  boundary  conditions  of  the  rise 
determine  the  position  of  the  first  and  the  last  points. 
Similarly,  velocity  boundary  conditions  at  both  ends  decide 
the  y-coordinates  of  the  second  and  the  second  last  control 
points.  Higher  derivative  boundary  conditions  will  have  a 
similar  effect  on  the  control  points  as  explained  in  Chapter 
II. 

From  the  above  discussion,  the  vector  of  design  varia- 
bles b,  can  be  expressed  as  follows: 

[b]  =  [Ro,  e,  bt,  PRlx,  PR2x PR9x,  PR4y,  PR5y,  PR6y, 

prix»  Pr2x,----,  Pr9x,  Pr4y,  Pr5y,  Pr6y  ]T 


(5.1) 


where 


R0  is  the  base  radius 

e  is  the  offset  of  roller  follower 

75 


bt  is  the  orientation  of  the  cam  blank  (This  design 
variable  is  applicable  only  for  the  minimization  of 
cutting  time) 

pRix  are  the  x-coordinates  of  the  control  points  for 

the  rise  cycle 
PRjy  are  the  y-coordinates  of  the  control  points  for 

the  rise  cycle 

prix  are  tne  x-coordinates  of  the  control  points  for 

the  return  cycle 
Prjy  are  the  y-coordinates  of  the  control  points  for 

the  return  cycle 


5.2   Cost 

The  minimization  of  a  cost  function  is  the  main  objec- 
tive of  optimization.  We  consider  two  different  classes  of 
minimization  problems: 

1)  Minimize  the  cutting  time:   The  objective  of  this 
minimization  problem  is  to  reduce  the  cutting  time 
required  to  manufacture  the  cam.   The  equations  used  to 
estimate  the  cutting  time  are  shown  in  Chapter  IV. 

2)  Minimize  the  cam  area:   In  many  cases,  the  material  cost 
is  a  more  important  factor  than  the  machining  time  or 
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machining  cost.  In  such  cases,  we  consider  the  volume 
of  material  required  to  be  the  cost  function.  As  dis- 
cussed in  Chapter  IV,  the  volume  of  material  required  is 
directly  proportional  to  the  cross-sectional  area  of  the 
cam.  Therefore,  the  cam  area  is  used  as  the  cost  func- 
tion. The  relevant  equations  for  calculation  of  cam 
area  are  shown  in  Chapter  IV. 


5 . 3  Constraints 

Defining  a  sufficient  set  of  constraints  is  an  impor- 
tant stage  in  optimal  design.  Constraints  provide  a  means 
of  preventing  the  optimization  process  from  converging  to  an 
unacceptable  solution.  Therefore,  the  designer  must  study 
the  system  carefully,  considering  all  performance  require- 
ments and  restrictions,  and  then  define  constraints  accord- 
ingly. 

In  this  thesis,  the  constraints  discussed  below  were 
imposed  as  standardized  constraints  for  all  problems, 
additional  problem-dependent  constraints  can  be  added  as 
needed. 
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1)  Elimination  of  follower  oscillation 

The  design  of  a  cam  is  generally  required  to  produce 
monotonic  rise  throughout  the  rise  cycle,  and  monotonic 
descent  throughout  the  return.  Failure  to  satisfy  these  re- 
quirements results  in  follower  oscillation  during  the  mo- 
tion. Oscillation  can  be  avoided  by  restricting  the  veloci- 
ty to  be  positive  throughout  the  rise  and  negative  through- 
out the  return.  These  requirements  can  be  enforced  through 
the  following  set  of  constraints: 

Rise  Cycle: 

g  =  -  Y'  <  0  (5.2) 

where 

Y'  is  the  reduced  velocity  of  the  follower 

Referring  to  Equations  (2.30),  (2.31)  and  (2.36),  Y'  can  be 
written  as 


n 
L  *  2  Piy  Ni/k'(u) 
i=0 


s  pix  Nj.,k'(u) 
i=0 


(5.3) 
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where 


L  is  the  total  lift 

pix  are  the  x-coordinates  of  the  control  points 
piy  are  the  y-coordinates  of  the  control  points 
R0  is  the  base  radius 
Ni,k(u)  are  the  blending  functions 

Ni,k'(u)  are  the  first  derivative  of  the  blending 
functions 


Return  Cycle: 


g  =  y>  <  o  (5-4) 


2)  Pressure  angle 

To  ensure  a  safe  design,  it  is  necessary  to  limit  the 
maximum  pressure  angle  to  be  within  a  specified  range.  This 
is  enforced  by  the  following  constraint  equation: 

g  =  02  -  0(J2  <  0  (5  5) 

where 

<p   is  the  instantaneous  pressure  angle. 

0d  is  the  specified  maximum  pressure  angle. 
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Substituting  from  Equation  (3.4),  the  expanded  form  of  the 
pressure  angle  constraint  is  given  by 

g  =  [  tan-1  (  (V  -  e)/(a+Y)  )  ]2  -  <pd2   <  0      (5-5) 

where 

e  is  the  offset  of  the  roller  follower 

a  is  the  vertical  distance  from  the  cam  axis  to  the 

base  circle,  given  by  Equation  (3.3) 
Y  is  the  lift,  given  by  Equation  (2.28) 
Y1  is  the  reduced  velocity,  given  by  Equation  (5.3) 

3 )  Undercutting 

The  curvature  of  the  cam  profile  must  be  adequate  to 
ensure  proper  follower  motion  and  to  prevent  excessive 
contact  stress  between  the  follower  and  the  cam.   This  is 
achieved  by  restricting  the  cam  curvature  to  be  greater  than 
one  and  a  half  times  of  the  radius  of  the  follower. 

g  =  1.5  *  rf  -  rc  <  o  (5.7) 

where 

rf  is  the  radius  of  the  roller  follower 
rc  is  the  instantaneous  radius  of  curvature  of  the  cam 
profile 
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Substituting  for  the  radius  of  curvature  of  the  profile,  we 
get 

[  C2  +   (Y'   -  e)2   ]3/2 
g  =  1.5*  Rf  -   <  o 

C2  +  (Y1  -  e) (2Y'  -  e)  -  C*Y' ' 

(5.8) 

where 

C  is  the  distance  between  the  cam  center  and  the 

profile,  given  by  Equation  (3.8) 
Y ' '  is  the  reduced  acceleration 

Referring  to  Equations  (2.30),  (2.31),  (2.32),  (2.33)  and 
(2.37),  the  reduced  acceleration  for  a  B-spline  is 


n                   n 
.s  pix  Ni/k'(u)  *  L*  S  Piy  Ni  k"(u) 
i=0                i=o 
Y«  = , 


.2  PixNi/k'(u) 

L  i=0  -I 

n  n 

L*   S   Piy  Ni  k'(u)  *  S  Pix  Ni#k"(u) 

1=0  i=o 

.S  pix  Ni,k'(u) 

L   1=0  J 
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(5.9) 


4)  Separation 

Separation  must  be  avoided  to  prevent  unacceptable  cam 
motions  as  well  as  other  problems  such  as  impact,  noise, 
etc.  This  can  be  achieved  by  restricting  the  driving  force 
to  be  positive  throughout  the  motion: 

9  =  "  Fd  ^  0  (5.10) 

where  F<j  is  the  driving  force,  given  by: 

Fd  =  (PL  +  Mf  Y)/  cos(0)  (5.11) 

where 

PL  is  the  specified  preload 

Mf  is  the  mass  of  the  roller  follower 

Y   is  the  true  acceleration 

The  true  acceleration  is  given  by 

5  "  Y"  *w2  (5.12) 

where  u  is  the  angular  velocity  of  cam. 

5)  Contact  Stress 

Excessive  stress  in  the  region  of  contact  between  the 
cam  and  the  follower  can  cause  surface  wear  and  failure  of 
the  cam  system.   This  can  be  avoided  by  restricting  the 
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g  =  Sc2/sd2  -  1.0  S   0  (5.13) 

where 

Sc  is  the  instantaneous  contact  stress,  given  by 
Equation  (3.6) 

Sd  is  the  specified  allowable  stress 

During  the  design  process,  it  is  important  to  ensure 
that  the  design  satisfies  all  the  specified  constraints. 
Each  rise  and  return  zone  is  divided  into  equal  intervals, 
and  the  checking  of  the  constraints  is  done  at  the  nodes 
between  these  intervals.  Considering  the  computing  time, 
the  distance  between  two  successive  nodes  must  not  be  too 
small,  otherwise  the  evaluation  of  constraint  functions  may 
be  too  expensive.  However,  if  the  interval  is  too  large,  it 
is  possible  to  miss  some  violations.  In  this  work,  each 
rise  zone  and  return  zone  is  divided  into  70  intervals.  In 
a  dwell  zone,  the  value  of  all  constraint  functions  is  con- 
stant. Therefore,  it  is  not  necessary  to  divide  the  dwell 
zones  into  intervals  to  check  for  constraint  violations. 
However,  the  behavior  of  constraints  in  a  dwell  zone  is 
accounted  for  by  evaluating  the  constraints  at  the  bound- 
aries of  the  rise  zone  and  the  return  zone. 
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For    other   design    problems,     additional    design    con- 
straints  may  be   needed.       In   computer   implementation   of   the 
method    presented    in    this    thesis,     additional    design    con- 
straints  can   be   added   conveniently   through   user   supplied 
codes. 

5.4     Optimization  Method 

After  the  optimal  design  problem  is  formulated,  there 
are  several  algorithms  available  to  solve  the  problem.  One 
class  of  solution  methods  is  called  sequential  unconstrained 
minimization  techniques  (SUMT) .  These  methods  find  a  solu- 
tion to  the  constrained  problem  by  solving  a  sequence  of 
unconstrained  problems.  The  sequence  of  solutions  obtained 
for  these  unconstrained  problems  converges  to  the  con- 
strained solution  [6].  This  class  of  methods  includes  the 
exterior  penalty   function  method  described  below. 

The  exterior  penalty  function  method  constructs  the 
unconstrained  problem  of  minimizing  the  following  "pseudo- 
objective   function"    : 

n(b,rp)    =  F(b)    +  rp  *  P(b)  (5.14) 

where 

F(b)  is  the  original  cost  function 
rp  is  the  penalty  multiplier 
P(b)  is  the  penalty  function 
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The  penalty  function,  P(b)  is  given  by 


m  1 

P(b)  =  S  {  max  [  0  ,  g-j  (b)  ]  }2  +  z   {  hk(b)  }2 
j=l  k=l 


(5.15) 


where 

gj (b)  are  the  inequality  constraints 
hfc(b)  are  the  equality  constraints 

From  Equations  (5.14)  and  (5.15),  we  observe  that  a  penalty 
is  imposed  if  and  only  if  one  or  more  constraints  are  vio- 
lated. When  violations  occur,  the  effect  of  the  penalty 
multiplier  is  to  magnify  these  violations  and  add  them  to 
the  true  cost,  thus  increasing  the  value  of  the  pseudo- 
objective  function.  Therefore,  the  unconstrained  optimiza- 
tion process  will  tend  to  move  away  from  the  infeasible 
regions,  since  constraint  violations  in  the  infeasible 
regions  cause  a  significant  increase  in  the  pseudo-objec- 
tive function  value. 

In  order  to  avoid  a  numerically  ill-conditioned  prob- 
lem, only  a  moderate  value  of  the  penalty  multiplier  is  used 
in  the  initial  optimization  stages,  and  the  value  is  in- 
creased as  the  optimization  progresses. 

For  the  minimization  of  cutting  time,  the  cutting  time 
is  a  discontinuous  function  of  the  design  variables.   This 
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is  because  the  cutting  time  is  dependent  on  the  number  of 
turns  of  the  tool  path.  A  small  change  in  the  value  of  a 
design  variable  may  cause  a  step  change  in  the  number  of 
turns.  Therefore,  a  gradient-based  method  is  not  applicable 
for  the  minimization  of  cutting  cost.  Consequently,  the 
unconstrained  minimization  method  selected  for  cutting  time 
minimization  was  a  directed  grid  search  [12].  This  is  a 
non-derivative  method  and  is  rather  crude,  but  it  performed 
satisfactorily  in  this  case. 

It  is  possible  to  implement  derivative-based  methods 
for  minimization  of  cam  area.  This  is  because  the  partial 
derivatives  of  cam  area  and  the  constraint  functions  with 
respect  to  the  design  variables  are  continuous  functions. 
The  Fletcher-Reeves  method  [6]  was  used  to  perform  the 
unconstrained  minimization  in  this  case.  This  method  re- 
quires the  first  order  partial  derivatives  of  the  pseudo- 
objective  function  with  respect  to  the  design  variables: 

tl'(b)  =  F'(b)  +  rp  *  P'(b)  (5.16) 

where 

fi' (b)  is  the  first  derivative  of  the  pseudo-objective 

function 
F' (b)  is  the  first  derivative  of  the  cost  function 
P'(b)  is  the  first  derivative  of  the  penalty  function 
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P1 (b)  is  given  by 

m 
P'b)  =  2*  E  max  [  0  ,  g-i(b)  ]  *  g-i  '  (b) 
1-1 

1 
+   E  2*  hk(b)  *  h'k(b)  (5.17) 

k=l 


where 

gj ' (b)  are  the  first  order  partial  derivatives  of 

inequality  constraints  g^  with  respect  to  the  design 

variables,  b 
hk' (b)  are  the  first  order  partial  derivatives  of  the 

equality  constraints  h^  with  respect  to  the  design 

design  variables,  b 

As  shown  in  Equation  (5.17),  in  order  to  use  Fletcher- 
Reeves  method,  it  is  necessary  to  evaluate  the  first  degree 
partial  derivatives  of  the  cost  and  constraint  functions. 
The  evaluation  of  these  derivatives  is  discussed  in  Chapter 
VI. 


87 


CHAPTER  VI 
DESIGN  SENSITIVITY  ANALYSIS 

In  chapter  V,  it  was  stated  that  a  gradient-based 
method  is  used  for  the  minimization  of  cam  area.  The 
Fletcher-Reeves  method  selected  for  this  work  requires  the 
evaluation  of  the  derivatives  of  cost  and  constraint  func- 
tions with  respect  to  the  design  vector  b  defined  in  Equa- 
tion (5.1)  . 

The  process  of  evaluating  these  derivatives  is  called 
"design  sensitivity  analysis"  and  is  discussed  in  detail  in 
this  chapter.  Only  the  derivation  of  non-zero  derivatives 
is  presented  here.  Any  derivative  that  does  not  appear  in 
the  following  discussion  is  a  derivative  with  a  constant 
value  of  zero. 


6.1  Design  Sensitivity  Analysis  of  Constraints 

1)  Elimination  of  oscillation 
Rise 

g  =  -  v   <  o 


<Sg  L    *    Ni/k'  (u)     *    D2 

7T in i=  i'9     t6-1' 

«P    y  Dl  « 
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£g_     l  *  Nj<k'(u) 

6Pjy         D1 


Where 


n 


n 


g  =  Y*  <  0 


«g       L  *  Ni|k'  (u)  *  D2 


4Pix  DX2 

«g      L  *  Nj/k'(u) 


*Pjy 


2)  Pressure  angle 


g  =  <t>2   -  0d2   <   0 


j    =4,6  (6.2) 


D2    =   2    Piy   Mi    k'(u)  (6.3) 

i=0 


Dl    =    S    Pix   Ni/k'(u)  (6.4) 

i=0 


Return 


i=    1,9  (6.5) 


j    =4,6  (6.1 
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From  equation  (5.7),  let 
C  =  a  +  Y 
V  =  (Y1  -  e)/C 


(6.7) 
(6.8) 


The  partial  derivatives  of  g  with  respect  to  design  varia- 
bles are  given  in  the  following  : 


Sq  2*0    SV 

SR0  (1+V2)    SRQ 

Sq  2*0        SV 


Se  (1+V2)       Se 

Sq  2*0       SV 


«Pix     (1+V2)     Spix 
Sq  2*0        SV 


SV 


3Y 


(1+V2)     «pj5 


1  =  1,9 


J  =  4,6 


(6.9) 


(6.10) 


(6.11) 


(6.12) 


Where 

SV 

= 

(e  -  Y' 
C2 

)    sc 

SR0 

SRQ 

SV 

= 

[C  + 

(Y'-e)] 

* 

C2 

<SC 

Se 

(Se 

SV 

SY> 

1 

•SPiv     5Pi 


(6.13) 


(6.14) 


(6.15) 
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SV  r  SY<  SC 

=    C  * (Y'  -  e)  *  

■5P-jv    L      <5P-jv  <SP 


C   * (Y'    -   e)    *   *   (6.16) 


<r3y  orjy  0Fjy 

SC  R„ 


5R0 

sqrt(R02  -  e2) 

SC 

e 

Ss 

sqrt(R02  -  e2) 

SC 

SPiY   " 

L  *  Nj/k(u) 

SY' 

L  *  »i,k' (U)  *  D2 

■SPix 

Dl2 

SY' 

L  *  Nj,k'(u) 

*pjy 

Dl 

(6.17) 


(6.18) 


(6.19) 


i=    1,9  (6.20) 


j    =4,6  (6.21) 


3 )  Undercutting 

g=1.5*rf-rc<0 

Referring  to  equation  (5.9),  let 

N  =  [  C2  +  (Y'  -  e)2  ]3/2  (6.22) 

D  =  C2  +  (Y'  -  e)*(2Y'  -  e)  -  C*Y"  (6.23) 
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sc  sc 

N    *     (2C    -    Y")*    -    3C    *    *[C2+(Y'-e)2]l/2*D 

*9  <5R0  SRa 


«R0  D2 

(6.24) 


r                                5C 
N    *  (2C    -    Y'1)     *    -    3*Y'    +       2*6 

«g  L  Se 


Ss  D2 


3D*[C2    +    (Y'-e)2]V2    *    |    c    * 


«C 

-    Y'    +    e 

Se 


D2 


SD  <SN 

N    *    -    D    *    

S<3      m  SPix SPix 

SPix  v2 


(6.25) 


i=    1,9  (6.26) 


Where 


=      3    [    C2    +    (Y1    -    e)2    ]-s    *    (Y'    -   e)     *   


SPi 


(6.27) 
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■  *  (4Y1  -  3e)   -  C  *  (6.28) 

spix        s?ix  6Pix 

Where 

*Y"  r  SB 


r  4N  1 

L   *        Di    *   -      3*N*   Niik'(u)         /(D!4) 


■SPix  L  tfp 


(6.29) 


T                                      <5C  5Y' 

N   *         (2C   -   Y")    *   +    (4Y"-3e)*   —      -   C* 

L  XT3  .  ct^  . 


if . ffjy *pjy  *PjyJ 

■SPjy  D2 


[5C  (SY ' 

C   *   +    (Y'    -   e    )* 

D2 


ipjy 

5PjyJ 

j=    4,6 

(6.30) 

Where 


5Y"    _       L    *    [    Dl    *    Nj,k"(u)     -    Ex    *    Nj(k'(u)     ] 
«Pjy  Dl3 


(6.31) 


93 


Where 


n 
El  =   2  Pix  Ni/k"(u)  (6.32) 

i=0 


4)  Separation 
g  =  -  Fd  <  0 

F(j  is  given  by  equation  (5.10). 

(   PL  +  MfY   )    *   sin(0   )    *  

49  5Rn 


*Ko  COS2(/> 

(    PL   +   MfY    )     *    sin(0    )     *   

S<3  Se 


COS20 


(6.33) 


(6.34) 


5<p 


MfUKi*    cos(0    j     *   +    (pL   +   MfY)*sin((/,    ,* 

Sq  &?ix  5Piv 


s?ix  cos2<p 


i=    1,9  (6.35) 


,                                  «"                                ..                            <S0 
"f1*   cos(0    )     *   +    (pL   +   MfY)*sin(0    )* 

£L      = ^JY iPjy 

*pjy  cos20 


j=    4,6  (6.36) 
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Referring   to   the    design    sensitivity    analysis    of   pressure 
angle  constraint,    we  have 


8<t>  SV  l 


<SR0  5R0  (1+V2) 


60  SV 


Se  Se  (1+V2) 


5<p  SV  X 


5pix  SPix  (1+V2) 

S<p  SV  1 

<SPjy  5Pjy  (1+V2) 

5)  Contact  Stress 


Si2 


Sd2 


-  1  <  0 


(6.37) 


(6.38) 


i  =1,9     (6.39) 


j  =4,6     (6.40) 


Referring  to  equation  (3.6),  let 

N  =  Fd  *  (rc  +  rf)  (6>41) 
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D  =  n    *  tc  *  rf  *  rc  *  [  (1-MC2)/EC  +  (1-Mf2)/Ef  ] 

(6.42) 

Dk  =  7T  *  tc   *  rc  *  [  (X-MC2)/EC  +  (1-Mf2)/Ef  ] 

(6.43) 


(6.44) 


Sq  1  5Si2 

■SRo    Sd2      5R0 

<sg  l  ssi2 

7T  =  ^*  "77  (6-45) 


4g  i  (SSi 


.2 


4pix      Sd2  SVi 


«g  l  «Si2 

<SPjy      Sd2  4Pjy 


i    =1,9  (6.46) 


j    =4,6  (6.47) 


Where 


*^c  SFd  -,  Sr, 

5S±2  SR0 


r                4rc                                      <SFd  -, 
D     Fd  *  +    (rc+rf)    *  -n*-      -*DV 

■SRo 


D2 

(6.48) 
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r  Src  SFd  -I  «rc 

,      -  D      Fd    *      +    (rc+rf)    *   -   N    *   *    Dk 

*Si2  L  <Se  se  J  <Se 


D2 


£«[a  4Fd  -j  <5rc 

isi2  *Pix 


(6.49) 


D     Fd   *     +    (rc+rf)    *  -  N  *  *   Dk 


D2 


i    =1,9  (6.50) 


r  Src  SFd  -I  Src 

.     -  D     Fd   *     +    (rc+rf)    *  -  n*  *   Dk 

!ii2    =         L  SPJY *FjyJ  *Pjy  k 


SPjy 


D2 


j    =4,6  (6.51) 

The  partial   derivatives   of  Fd,    with  respect  to  the  design 
variables,     can    be    obtained    from    the    design    sensitivity 
analysis  of   separation  constraint. 


(    PL  +   MfY    )    *    sin(0    )    *   

^d  <SR0 

,„      *  : (6.52) 

oRq  cos'0 
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(  PL  +  MfY  )  *  sin(0  )  *  

"   ' (6.53) 

Se  cos20 


Mf  w**   cos(0  )  *  +  (pL  +  MfY)*sin(0  )* 

SFd SPix <5Pj 

<SPix  cos20 


i=  1,9  (6-54) 


Mf  u^*  cos(0  )  *  +  (pL  +  MfY)*sin(0  )* 

£1*    . 5l21 ffjy 

<5Pjy  COS20 


j*  4,6  (6.55) 

From  the  design  sensitivity  analysis  of  undercutting  con- 
straint, we  have 


SC 

N  *  (2C  -  Y1 ' ) *  


5R0  D2 


<5C 

3C  *  *[c2+(Y'-e)2]V2*D 

5P-o 

D2 
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(6.56) 


5rr 


Se 


N   * 


(2C   -   Y' 


sc 

Se 


3*Y"    +      2*e 


D2 


3D*[C2   +    (Y'-e)2]l/2    *    \    c 


D2 


<5C 
*   -   Yp    +   e 


(6.57) 


Srr 


SPi 


<5D 


5N 


N    * 


D    * 


<SPi 


*Pi 


V2 


i=    1,9  (6.58) 


«p 


r  <SC  <sy'  8t< 

N*|^(2C    -    Y")     *   -     -   +    (4Y"-3e)*   -         -    C*- 


<SP 


JY 


<SP 


3Y 


5P 


:y 


3y 


D2 


3D*[C2    +    (Y'-e)2]V2    *        c   *   -    -  +    (Y'    -   e    ) *   


5P 


]y 


<5P 


]y 


D2 


j=    4,6  (6.59) 
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6.2   Design  Sensitivity  Analysis  of  Cost  Function 
Cam  Area 

The  equations  to  evaluate  cam  area  were  given  in  Chap- 
ter IV.  From  those  equations,  the  partial  derivatives  of 
the  cam  area  with  respect  to  the  design  variables  can  be 
derived  as  follows: 


Sk  fui 


(R0  +  L*Y)  *  Di 
•SRo    .Wo 


du 


+  R0*eb  +  (R0+L)*et 


(6.60) 


5A 


«Pi 


=   0.5  * 


(R0  +  L*Y)2   *  Nifk' (u)  du 


i  =1,9 


(6.61) 


<5A 


.[»' 


■SPjy     JUC 


L  *  (RQ  +  L*Y)   *  Njfk(u)*  D1   du 


j  =4,6 


(6.62) 
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CHAPTER  VII 
NUMERICAL  EXAMPLES 

The  technique  for  cam  design  using  B-spline  curves  as 
described  in  the  preceding  chapters  was  implemented  in  a 
computer  program.  Several  numerical  examples  were  solved  to 
test  the  computer  program  and  to  verify  the  design  technique 
developed  in  this  work. 

The  physical  properties  for  the  material  of  cam  and  the 
follower  were  taken  from  [13].  These  properties  were  kept 
unchanged  in  all  the  examples. 

Physical  Properties 

Material:  Hardened  Steel  SAE  1010  (for  both  cam  &  follower) 
Density  of  the  material  =  0.28  lb/in2 
Young's  modulus  of  elasticity  =  29  x  106  psi 
Poison's  ratio  =  0.29 

7.1  Minimization  of  cam  cutting  Time 

The  machining  data  required   for  this  problem  was 
obtained  from  [10].   These  parameters  were  also  kept  un- 
changed in  all  the  examples  presented  in  this  chapter. 
Machining  Data: 

Tool  Material:  High  Speed  Steel 

Tool  diameter  =  0.75  inches 

Number  of  teeth  =  4 

Cutting  speed  =  109.95  fpm 

Feed  per  tooth  =  0.0025  inch/tooth 

Spindle  speed  =  560  RPM 

Feed  rate  =  9.333  x  10"2  inch/sec 
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The  first  four  examples  presented  here  treat  the  cut- 
ting time  as  the  cost  function  to  be  minimized.  The  design 
variables  used  in  the  cutting  time  minimization  problems  are 
given  below: 

Number  of  design  variables  =  27 

B(l)  :  the  base  radius  of  the  cam 

B(2)  :  the  offset  of  the  roller  follower 

B(3)  :  the  orientation  of  the  blank  on  the  NC  table 

B(4)  :  the  x-coordinate  of  the  second  control  point  for 

the  rise  cycle. 
B(5)  :  the  x-coordinate  of  the  third  control  point  for 

the  rise  cycle. 
B(6)  :  the  x-coordinate  of  the  fourth  control  point  for 

the  rise  cycle. 
B(7)  :  the  x-coordinate  of  the  fifth  control  point  for 

the  rise  cycle. 
B(8)  :  the  y-coordinate  of  the  fifth  control  point  for 

the  rise  cycle. 
B(9)  :  the  x-coordinate  of  the  sixth  control  point  for 

the  rise  cycle. 
B(10)  :  the  y-coordinate  of  the  sixth  control  point  for 

the  rise  cycle. 
B(ll)  :  the  x-coordinate  of  the  seventh  control  point  for 

the  rise  cycle. 
B(12)  :  the  y-coordinate  of  the  seventh  control  point  for 

the  rise  cycle. 
B(13)  :  the  x-coordinate  of  the  eighth  control  point  for 

the  rise  cycle. 
B(14)  :  the  x-coordinate  of  the  ninth  control  point  for 

the  rise  cycle. 
B(15)  :  the  x-coordinate  of  the  tenth  control  point  for 

the  return  cycle. 
B(X6)  :  the  x-coordinate  of  the  second  control  point  for 

the  return  cycle. 
B(17)  :  the  x-coordinate  of  the  third  control  point  for 

the  return  cycle. 
B(18)  :  the  x-coordinate  of  the  fourth  control  point  for 

the  return  cycle. 
B(19)  :  the  x-coordinate  of  the  fifth  control  point  for 

the  return  cycle. 
B(20)  :  the  y-coordinate  of  the  fifth  control  point  for 

the  return  cycle. 
B(21)  :  the  x-coordinate  of  the  sixth  control  point  for 
the  return  cycle. 
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B(22)  :  the  y-coordinate  of  the  sixth  control  point  for 

the  return  cycle. 
B(23)  :  the  x-coordinate  of  the  seventh  control  point  for 

the  return  cycle. 
B(24)  :  the  y-coordinate  of  the  seventh  control  point  for 

the  return  cycle. 
B(25)  :  the  x-coordinate  of  the  eighth  control  point  for 

the  return  cycle. 
B(26)  :  the  x-coordinate  of  the  ninth  control  point  for 

the  return  cycle. 
B(27)  :  the  x-coordinate  of  the  tenth  control  point  for 

the  return  cycle. 
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Example  1; 

Preload  =  320.0  lbs 

Maximum  allowable  stress  =  250,000  psi 

Maximum  allowable  pressure  angle  =25.0  degrees 

Angular  velocity  of  cam  =  800  RPM 

Mass  of  follower  mechanism  =  0.02  slugs 

Lift  =  1.0  inches 

Cam  thickness  =  0.25  inches 

Follower  diameter  =  0.85  inches 

Follower  thickness=  0.25  inches 

Required  accuracy  on  the  profile  =  0.01  inches 


Cycle  Characteristics: 

Rise         Dwell  Return 

0.0-170.0      170.0-190.0      190.0-340.0 


Dwell 
340.0-360.0 


Design  Variable 

B(l) 

B(2) 

B(3) 

B(4) 

B(5) 

B(6) 

B(7) 

B(8) 

B(9) 

B(10) 

B(ll) 

B(12) 

B(13) 

B(14) 

B(15) 

B(16) 

B(17) 

B(18) 

B(19) 

B(20) 

B(21) 

B(22) 

B(23) 

B(24) 

B(25) 

B(26) 

B(27) 


Initial  Value 


4 

.0000 

0 

.0000 

0 

.2094 

0 

.4273 

0 

.4987 

0 

.5700 

1 

.3532 

0 

.1300 

1 

.7805 

0 

.6500 

2 

.2655 

0. 

.9000 

2 

.3460 

2 

.4956 

2. 

.9308 

3. 

,8729 

3. 

.9211 

3. 

9694 

4. 

4970 

1. 

0000 

5. 

0253 

0. 

3000 

5. 

2653 

0. 

1000 

5. 

5053 

5. 

7453 

5. 

8653 

Final  Value 

3 

.6400 

0 

.0500 

0 

.7540 

0 

.3827 

0 

.4245 

0 

.8668 

1 

.0564 

0 

.1000 

1 

.4837 

0 

.6700 

1 

.8501 

1 

.0000 

2, 

.2125 

2. 

,3621 

2. 

,4561 

3. 

,3230 

3. 

3583 

3. 

4458 

4. 

4184 

0. 

7700 

5. 

0253 

0. 

4700 

5. 

2129 

0. 

1400 

5. 

4791 

5. 

5489 

5. 

7475 
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Comment 

The  design  procedure  started  within  the  infeasible 
region  with  the  maximum  violation  of  484.858  for  the  separa- 
tion constraint.  When  the  design  stepped  into  the  feasible 
region  for  the  first  time,  the  base  radius  and  cutting  time 
were  found  to  be  5.80  inches  and  1607.48  seconds  respective- 
ly. Note  that  the  jerk  for  the  return  cycle  is  actually 
started  at  zero  as  shown  in  Figure  7.2. 


Maximum  Values 

Force 

Pressure  Angle 
Contact  Stress 
Reduced  Velocity 
Reduced  Acceleration 
Reduced  Jerk 


Magnitude 

459.051  lbs 

-16.863' 

1.708E+5  psi 
1.196  in/s 
3.914  in/s2 
-212.666  in/s3 


Location  On  Cam  Profile 
(degree) 

58.053 
295.431 

58.053 

73.194 

56.388 
190.194 


Highest  constraint  violation  =  1.480  E-6 

Total  number  of  constraint  violated  =  1 

Number  of  turns  =  1530 

Tool  path  distance  =  3  3.00  inches 

Cutting  time  =  965.55  seconds 

%  reduction  in  cutting  time  from 
the  first  feasible  design  =  39.93 
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RISE    CYCLE 


o      3.500 


(X10  °°) 


RISE    CYCLE 


(XI  o  °°) 


Figure   7 . 1 


Performance  Characteristics 
for  Example  1  (Rise  Cycle) 
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RETURN    CYCLE 


(X10   °°) 


RETURN    CYCLE 


(XI  o  °°) 


Figure  7.2  :  Performance  Characteristics 
for  Example  1  (Return  Cycle) 
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. . .  Cam  profile  of  initial  design 

Cam  profile  of  final  design 

Base  circle  of  final  design 


Figure  7.3  :  Cam  Profiles  for  Example  1 
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Example  2: 

Preload  =  370.0  lbs 

Maximum  allowable  stress  =  250,000  psi 

Maximum  allowable  pressure  angle  =  25.0  degrees 

Angular  velocity  of  cam  =  600  RPM 

Mass  of  follower  mechanism  =  0.02  slugs 

Lift  =  1.3  inches 

Cam  thickness  =  0.2  5  inches 

Follower  diameter  =  0.85  inches 

Follower  thickness=  0.25  inches 

Required  accuracy  on  the  profile  =0.01  inches 


Cycle  Characteristics : 

Rise  Dwell  Return 

0.0-100.0      100.0-160.0      160.0-320.0 


Dwell 
320.0-360.0 


Design  Variable 

B(l) 

B(2) 

B(3) 

B(4) 

B(5) 

B(6) 

B(7) 

B(8) 

B(9) 

B(10) 

B(ll) 

B(12) 

B(13) 

B(14) 

B(15) 

B(16) 

B(17) 

B(18) 

B(19) 

B(20) 

B(21) 

B(22) 

B(23) 

B(24) 

B(25) 

B(26) 

B(27) 


Initial  Value 


3 

.5000 

0 

.0000 

0 

.2094 

0 

.2513 

0 

.2933 

0 

.3353 

0 

.7960 

0 

.1300 

1 

.0473 

0 

.6500 

1 

.3327 

0 

.9000 

1. 

.3800 

1 

.4680 

1. 

.7240 

3. 

,3863 

3. 

,4379 

3. 

,4894 

4. 

0521 

1. 

0000 

4. 

6156 

0. 

3000 

4. 

8716 

0. 

1000 

5. 

1760 

5. 

3836 

5. 

5116 

Final  Value 

5.2100 
0.1000 
0.0628 
0.1641 
0.2933 

.3528 

.6215 

.1200 

.8728 

.6600 

.1407 

.9100 

.3887 
1.5465 
1.6455 
3.1210 
3.1447 
3.7687 
4.0800 
0.7000 
4.5877 
0.0900 
4.8437 
0.0000 
5.1276 
5.3836 
5.5116 


0. 
0. 
0. 
0. 
0. 

1. 

0. 

1. 
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Comment 

The  design  procedure  started  within  the  infeasible 
region  with  the  maximum  violation  of  1722.565  for  the  under- 
cutting constraint.  The  separation  constraint  was  also 
violated  with  a  maximum  violation  of  1354.038.  When  the 
design  stepped  into  the  feasible  region  for  the  first  time, 
the  base  radius  and  cutting  time  were  found  to  be  6.2  0 
inches  and  1699.40  seconds  respectively. 


Maximum  Values        Magnitude      Location  On  Cam  Profile 

(degree) 

Force  497.570  lbs  33.972 

Pressure  Angle  19.557°  42.806 

Contact  Stress  1.541E+5  psi  36.908 

Reduced  Velocity        2.148  in/s  44.278 

Reduced  Acceleration    5.712  in/s2  32.520 

Reduced  Jerk  -58.623  in/s3  42.806 


Highest  constraint  Violation  =  0.000 

Total  number  of  constraint  violated  =  0 

Number  of  turns  =  2192 

Tool  path  distance  =  47.0  inches 

Cutting  time  =  1381.24  seconds 

%  reduction  in  cutting  time  from 
the  first  feasible  design  =  18.72 
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Figure  7.4  :  Performance  Characteristics 
for  Example  2  (Rise  Cycle) 
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RETURN  CYCLE 
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Figure  7.5  :  Performance  Characteristics 
for  Example  2  (Return  Cycle) 
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Cam  profile  of  initial  design 

Cam  profile  of  final  design 

-  -   Base  circle  of  final  design 


Figure  7.6  :  Cam  Profiles  for  Example  2 
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Example  3: 

Preload  =  380.0  lbs 

Maximum  allowable  stress  =  250,000  psi 

Maximum  allowable  pressure  angle  =28.0  degrees 

Angular  velocity  of  cam  =  800  RPM 

Mass  of  follower  mechanism  =0.02  slugs 

Lift  =  1.2  inches 

Cam  thickness  =  0.25  inches 

Follower  diameter  =  0.85  inches 

Follower  thickness=  0.25  inches 

Required  accuracy  on  the  profile  =0.01  inches 


Cycle  Characteristics : 

Rise  Dwell  Return 

0.0-150.0      150.0-180.0      180.0-330.0 


Dwell 
330.0-360.0 


Design  Variable 

B(l) 
B(2) 
B(3) 

B(4) 

B(5) 

B(6) 

B(7) 

B(8) 

B(9) 

B(10) 

B(ll) 

B(12) 

B(13) 

B(14) 

B(15) 

B(16) 

B(17) 

B(18) 

B(19) 

B(20) 

B(21) 

B(22) 

B(23) 

8(24) 

B(25) 

B(26) 

B(27) 


Initial  Value 


5 

.4000 

0 

.0000 

0 

.3491 

0 

.3770 

0 

.4400 

0 

.5030 

0 

.9322 

0 

.1300 

1 

.3092 

0 

.8500 

1 

.9990 

1. 

.0000 

2 

.0700 

2. 

.3333 

2. 

.4551 

3, 

5674 

3. 

6151 

3. 

7949 

4. 

0606 

0. 

9000 

4. 

8507 

0. 

2000 

5. 

0907 

0. 

1000 

5. 

3307 

5. 

5707 

5. 

6908 

inal  Value 

5 

.2200 

0 

.0500 

0 

.9283 

0 

.2592 

0 

.2960 

0 

.5292 

0 

.9322 

0 

.0000 

1 

.3092 

0 

.8500 

1 

.9990 

1 

.0000 

2 

.2140 

2 

.3400 

2. 

,4551 

3. 

.1878 

3. 

.2099 

3. 

.8211 

4. 

3486 

0. 

9900 

4. 

8769 

0. 

2900 

5. 

0907 

0. 

0000 

5. 

3308 

5. 

4268 

5. 

5599 
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Comment 


The  design  procedure  started  within  the  feasible  region 
with  an  initial  cutting  time  of  1531.54  seconds. 


Maximum  Values 

Force 

Pressure  Angle 
Contact  Stress 
Reduced  Velocity 
Reduced  Acceleration 
Reduced  Jerk 


Magnitude 

582.177  lbs 

-15.798° 

1.664E+5  psi 
1.645  in/s 
5.534  in/s2 
34.261  in/s3 


Location  On  Cam  Profile 
(degree) 

290.586 

277.835 

290.586 

62.450 

291.920 

284.866 


Highest  constraint  violation  =  2.004E-6 

Total  number  of  constraint  violated  ■  1 

Number  of  turn  =  2196 

Cutting  time  =  1377.60  seconds 

Tool  path  distance  =  46.59  inches 

%  Reduction  in  cutting  time  =  10.05 
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Figure  7.7  :  Performance  Characteristics 
for  Example  3  (Rise  Cycle) 
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Figure  7.8  :  Performance  Characteristics 
for  Example  3  (Return  Cycle) 
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Cam  profile  of  initial  design 
Cam  profile  of  final  design 
Base  circle  of  final  design 


Figure  7.9  :  cam  Profiles  for  Example  3 
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Example  4: 

Preload  =  400.0  lbs 

Maximum  allowable  stress  =  250,000  psi 

Maximum  allowable  pressure  angle  =  20.0  degrees 

Angular  velocity  of  cam  =  1000  RPM 

Mass  of  follower  mechanism  =  0.02  slugs 

Lift  =  1.0  inches 

Cam  thickness  =  0.25  inches 

Follower  diameter  =  0.85  inches 

Follower  thickness=  0.25  inches 

Required  accuracy  on  the  profile  =0.01  inches 


Cycle  Characteristics : 

Rise  Dwell  Return 

-20.0-130.0      130.0-140.0      140.0-310.0 


Dwell 
310.0-340.0 


Design  Variable 

B(l) 

B(2) 

B(3) 

B(4) 

B(5) 

B(6) 

B(7) 

B(8) 

B(9) 

B(10) 

B(ll) 

B(12) 

B(13) 

B(14) 

B(15) 

8(16) 

B(17) 

B(18) 

B(19) 

B(20) 

B(21) 

B(22) 

B(23) 

B(24) 

B(25) 

B(26) 

B(27) 


Initial  Value 


3 

.8000 

0 

.0000 

1 

.2566 

0 

.0279 

0 

.0909 

0 

.4157 

0 

.5831 

0 

.1300 

1 

.2219 

0 

.9500 

1 

.6499 

1 

.0000 

1 

.7209 

1, 

.8529 

1 

.9751 

2, 

,7776 

2. 

,8324 

2, 

8872 

3. 

7817 

0. 

9000 

4. 

3805 

0. 

3000 

4. 

6525 

0. 

0000 

4. 

9245 

5. 

1965 

5. 

3325 

in. 

al  Value 

2 

.4500 

0 

.0000 

1 

.3194 

0 

.0279 

0 

.1171 

0 

.2063 

0 

.6093 

0 

.1300 

1 

.2219 

0 

.9500 

1 

.7808 

0 

.9600 

1, 

.8126 

2, 

.0754 

2. 

,1846 

2. 

,7776 

2, 

8027 

2. 

8872 

3. 

4850 

0. 

8900 

4. 

3509 

0. 

4800 

4. 

6229 

0. 

0700 

4. 

9245 

4. 

9444 

5. 

0803 
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Comment 

The  design  procedure  started  within  the  infeasible 
region  with  violations  in  undercutting,  oscillation  and 
stress  constraints.  The  undercutting  constraint  was  most 
violated  with  a  maximum  violation  of  568.26.  When  the 
design  stepped  into  the  feasible  region  for  the  first  time, 
the  base  radius  and  the  cutting  time  were  found  to  be  equal 
to  4.7  inches  and  1257.63  seconds  respectively. 


Maximum  Values 

Force 

Pressure  Angle 
Contact  Stress 
Reduced  Velocity 
Reduced  Acceleration 
Reduced  Jerk 


Magnitude 

564.056  lbs 

-19.521° 

1.653E+5  psi 
1.016  in/s 
2.720  in/s2 

-44.969  in/s3 


Location  On  Cam  Profile 
(degree) 

269.931 

259.760 

266.744 

49.600 

269.930 

160.132 


Highest  constraint  violation  =  0.000000 

Total  number  of  constraint  violated  =  0 

Number  of  turn  =  1062 

Cutting  time  =  676.57  seconds 

Tool  path  distance  =  2  3.49  inches 

%   reduction  in  cutting  time  from 

the  first  feasible  design  =  46.20 
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Figure  7.10  :  Performance  Characteristics 
for  Example  4  (Rise  Cycle) 
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RETURN    CYCLE 
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Figure  7.11  :  Performance  Characteristics 
for  Example  4  (Return  cycle) 
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Cam  profile  of  initial  design 
Cam  profile  of  final  design 


Base  circle  of  final  design 


Figure  7.12  :  Cam  Profiles  for  Example  4 
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7.2  Minimization  of  Cam  Material 

The  next  three  examples  treat  the  cross-sectional  area 
of  the  cam  as  the  cost  function  to  be  minimized.  The  design 
variables  for  the  minimization  of  cam  material  are  given 
below: 

Number  of  design  variables  =26 

B(l)  :  the  base  radius  of  the  cam 

B(2)  :  the  offset  of  the  roller  follower 

B(3)  :  the  x-coordinate  of  the  second  control  point  for 

the  rise  cycle. 
B(4)  :  the  x-coordinate  of  the  third  control  point  for 

the  rise  cycle. 
B(5)  :  the  x-coordinate  of  the  fourth  control  point  for 

the  rise  cycle. 
B(6)  :  the  x-coordinate  of  the  fifth  control  point  for 

the  rise  cycle. 
B(7)  :  the  y-coordinate  of  the  fifth  control  point  for 

the  rise  cycle. 
B(8)  :  the  x-coordinate  of  the  sixth  control  point  for 

the  rise  cycle. 
B(9)  :  the  y-coordinate  of  the  sixth  control  point  for 

the  rise  cycle. 
B(10)  :  the  x-coordinate  of  the  seventh  control  point  for 

the  rise  cycle. 
B(ll)  :  the  y-coordinate  of  the  seventh  control  point  for 

the  rise  cycle. 
B(12)  :  the  x-coordinate  of  the  eighth  control  point  for 

the  rise  cycle. 
B(13)  :  the  x-coordinate  of  the  ninth  control  point  for 

the  rise  cycle. 
B(14)  :  the  x-coordinate  of  the  tenth  control  point  for 

the  return  cycle. 
B(15)  :  the  x-coordinate  of  the  second  control  point  for 

the  return  cycle. 
B(16)  :  the  x-coordinate  of  the  third  control  point  for 

the  return  cycle. 
B(17)  :  the  x-coordinate  of  the  fourth  control  point  for 

the  return  cycle. 
B(18)  :  the  x-coordinate  of  the  fifth  control  point  for 

the  return  cycle. 
B(19)  :  the  y-coordinate  of  the  fifth  control  point  for 
the  return  cycle. 
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B(20)  :  the  x-coordinate  of  the  sixth  control  point  for 

the  return  cycle. 
B(21)  :  the  y-coordinate  of  the  sixth  control  point  for 

the  return  cycle. 
B(22)  :  the  x-coordinate  of  the  seventh  control  point  for 

the  return  cycle. 
B(23)  :  the  y-coordinate  of  the  seventh  control  point  for 

the  return  cycle. 
B(24)  :  the  x-coordinate  of  the  eighth  control  point  for 

the  return  cycle. 
B(25)  :  the  x-coordinate  of  the  ninth  control  point  for 

the  return  cycle. 
B(26)  :  the  x-coordinate  of  the  tenth  control  point  for 

the  return  cycle. 


125 


Example  5 : 

Preload  =  340.0  lbs 

Maximum  allowable  stress  =  250,000  psi 

Maximum  allowable  pressure  angle  =18.0  degrees 

Angular  velocity  of  cam  =  800  RPM 

Mass  of  follower  mechanism  =  0.02  slugs 

Lift  =1.00  inches 

Cam  thickness  =  0.25  inches 

Follower  diameter  =  0.85  inches 

Follower  thickness=  0.25  inches 


Cycle  Characteristics 


Rise 

Dwell 

Return 

Dwell 

0.0-150.0      150 

.0-180.0      180.0-330. 

.0      330.0-360.0 

Design  Variable 

Initial  Value 

Final  Value 

8(1) 

4.5000 

3.7344 

B(2) 

0.0000 

0.0001 

B(3) 

0.3770 

0.3772 

B(4) 

0.4400 

0.4424 

B(5) 

0.5030 

0.5148 

B(6) 

1.1940 

1.2282 

B(7) 

0.1300 

0.0772 

B(8) 

1.5710 

1.5945 

B(9) 

0.6500 

0.6242 

B(10) 

1.9990 

1.9231 

B(ll) 

0.9000 

0.9094 

B(12) 

2.0700 

2.1345 

B(13) 

2.2020 

2.2648 

B(14) 

2.5860 

2.5848 

B(15) 

3.2860 

3.2859 

B(16) 

3.3440 

3.3433 

B(17) 

3.4020 

3.3922 

B(18) 

4.0350 

3.9976 

B(19) 

1.0000 

0.9390 

B(20) 

4.4500 

4.4062 

B(21) 

0.3000 

0.2445 

B(22) 

4.9570 

4.9349 

B(23) 

0.1000 

0.0544 

B(24) 

5.2450 

5.2382 

B(25) 

5.5330 

5.5314 

B(26) 

5.6770 

5.6768 
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Comment 

The  design  procedure  started  within  the  feasible  region 
with  an  initial  area  of  78.364  inches. 


Maximum  Values  Magnitude      Location  On  Cam  Profile 

(degree) 

Force  632.654  lbs           74.918 

Pressure  Angle  15.121'             88.806 

Contact  Stress  1.782E+5  psi        74.918 

Reduced  Velocity  1.153  in/s          88.806 

Reduced  Acceleration  -2.426  in/s2       120.756 

Reduced  Jerk  -20.991  in/s3       190.331 


Highest  constraint  violation  =  0.000000000 
Total  number  of  constraint  violated  =  0 
Area  =  55.493  inches 
%  reduction  in  area  =  29.186 
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Figure   7.13 


Performance  Characteristics 
for  Example  5  (Rise  Cycle) 
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Figure  7.14  :  Performance  Characteristics 
for  Example  S  (Return  Cycle) 
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Cam  profile  of  initial  design 
Cam  profile  of  final  design 
Base  circle  of  final  design 


Figure  7.15  :  Cam  Profiles  for  Example  5 
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Example  6: 

Preload  =  3  60.0  lbs 

Maximum  allowable  stress  =  250,000  psi 

Maximum  allowable  pressure  angle  =  2  5.0  degrees 

Angular  velocity  of  cam  =  650  RPM 

Mass  of  follower  mechanism  =  0.02  slugs 

Lift  =  1.50  inches 

Cam  thickness  =  0.25  inches 

Follower  diameter  =  0.85  inches 

Follower  thickness=  0.25  inches 


Cycle  Characteristics : 

Rise  Dwell  Return 

180.0-310.0    310.0-340.0      -20.0-110.0 


Dwell 
110.0-180.0 


Design  Variable 

B(D 

B(2) 

B(3) 

B(4) 

B(5) 

B(6) 

B(7) 

B(8) 

B(9) 

B(10) 

B(ll) 

B(12) 

B(13) 

B(14) 

B(15) 

B(16) 

B(17) 

B(18) 

B(19) 

B(20) 

B(21) 

B(22) 

B(23) 

B(24) 

B(25) 

B(26) 


Initial  Value 


4 

.2000 

0 

.0512 

3 

.4683 

3 

.5229 

3 

.5775 

4 

.1764 

0 

.1300 

4 

.4851 

0 

.6675 

4 

.7740 

0 

.9411 

4 

.9933 

5. 

.1116 

5 

.3814 

0, 

.1334 

0. 

,1753 

0. 

.2172 

0. 

6743 

1, 

0000 

1. 

1322 

0. 

3000 

1. 

3402 

0. 

1000 

1. 

5482 

1. 

7562 

1. 

8602 

Final  Value 

3.7500 

0.0516 
3.4685 

5248 

5872 

2050 

0920 

5204 
0.6337 
4.7953 
0.9120 
4.9998 
5.1127 
5.3818 

1334 

1746 

2084 

6416 

9560 

0947 

2653 

3218 
0.0769 
1.5427 
1.7549 
1.8600 
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Comment 

The  design  procedure  started  within  the  infeasible 
region  with  an  initial  area  of  81.8  68  inches.  The  separa- 
tion constraint  was  violated  with  a  maximum  violation  of 
824.50.   The  oscillation  constraint  was  also  violated. 

Maximum  Values        Magnitude      Location  On  Cam  Profile 

(degree) 

Force  824.503  lbs  71.578 

Pressure  Angle          23.183°  253.083 

Contact  Stress           1.992E+5  psi  68.673 

Reduced  Velocity         2.002  in/s  256.732 

Reduced  Acceleration     4.373  in/s2  71.578 

Reduced  Jerk  -86.491  in/s3  9.505 


Highest  constraint  violation  =  4.307E-6 
Total  number  of  constraint  violated  =  1 
Area  =  67.373  inches 
%  reduction  in  area  =  17.707 
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Figure  7.16  :  Performance  Characteristics 
for  Example  6  (Rise  Cycle) 
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Figure  7.17  :  Performance  Characteristics 
for  Example  6  (Return  Cycle) 
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Cam  profile  of  initial  design 
Cam  profile  of  final  design 


-   Base  circle  of  final  design 


Figure  7.18  :  cam  Profiles  for  Example  6 
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Example  7: 

Preload  =  3  60.0  lbs 

Maximum  allowable  stress  =  250,000  psi 

Maximum  allowable  pressure  angle  =20.0  degrees 

Angular  elocity  of  cam  =7  00  RPM 

Mass  of  follower  mechanism  =0.02  slugs 

Lift  =  l.io  inches 

Cam  thickness  =  0.25  inches 

Follower  diameter  =0.85  inches 

Follower  thickness=  0.2  5  inches 


Cycle  Characteristics : 
Rise  Dwell 

0.0-140.0      140.0-150.0 


Return 
150.0-330.0 


Dwell 
330.0-360.0 


Design  Variable 

B(l) 

B(2) 

B(3) 

B(4) 

8(5) 

B(6) 

B(7) 

B(8) 

B(9) 

B(10) 

B(ll) 

B(12) 

B(13) 

B(14) 

B(15) 

B(16) 

B(17) 

B(18) 

B(19) 

B(20) 

B(21) 

B(22) 

B(23) 

B(24) 

B(25) 

B(26) 


Initial  Value 


5 

.0000 

0 

.0000 

0 

.3519 

0 

.4107 

0 

.4695 

1 

.1144 

0 

.1300 

1 

.4484 

0 

.6687 

1 

.7657 

0 

.9442 

1, 

.9896 

2 

.1169 

2 

,4122 

3, 

,2860 

3. 

,3440 

3. 

4020 

4. 

0350 

1. 

0000 

4. 

6690 

0. 

3000 

4. 

9570 

0. 

1000 

5. 

2450 

5. 

5330 

5. 

6770 

Final  Value 

3 

.5655 

0 

.0001 

0 

.3523 

0 

.4148 

0 

.4918 

1 

.1822 

0 

.0301 

1 

.5272 

0 

.5940 

1 

.7711 

0 

.8981 

2, 

,0329 

2, 

.1459 

2, 

,4097 

3, 

.2858 

3, 

3415 

3. 

3794 

3. 

9570 

0. 

8614 

4. 

5813 

0. 

1865 

4. 

9146 

0. 

0208 

5. 

2327 

5. 

5302 

5. 

6767 

136 


Comment 

The  design  procedure  started  within  the  feasible  region 
with  an  initial  area  of  99.742  inches. 


Maximum  Values 


Force 

Pressure  Angle 
Contact  Stress 
Reduced  Velocity 
Reduced  Acceleration 
Reduced  Jerk 


Magnitude 


699.396  lbs 
20.153° 
1.861E+5  psi 
1.525  in/s 
-3.324  in/s2 
-93.409  in/s3 


Location  On  Cam  Profile 
(degree) 

76.552 

84.511 

76.552 

86.277 
191.603 
190.698 


Highest  constraint  violation  =  2.134E-5 
Total  number  of  constraint  violated  =  2 
Area  =  54.017  inches 
%  reduction  in  area  =  45.84  0 
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Figure  7.19 


Performance  Characteristics 
for  Example  7  (Rise  Cycle) 
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RETURN     CYCLE 
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Figure  7.20  :  Performance  Characteristics 
for  Example  7  (Return  Cycle) 
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Cam  profile  of  initial  design 
Cam  profile  of  final  design 
Base  circle  of  final  design 


Figure  7.21  :  cam  Profiles  for  Example 
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CHAPTER  VIII 
CONCLUSIONS 


In  this  thesis,  a  new  technique  for  designing  plate 
cams  has  been  developed.  The  approach  adopted  was  to  model 
the  cam  profile  by  B-spline  curves,  and  use  the  B-spline 
representation  in  the  design  process.  In  addition  to  satis- 
fying the  boundary  conditions,  B-splines  are  capable  of 
providing  design  control  over  the  shape  of  the  profile  in 
the  interior  of  a  motion  cycle,  which  is  a  crucial  factor 
for  satisfying  the  specified  design  requirements  in  this 
region.  B-splines  possess  several  attractive  properties  and 
can  be  used  to  model  a  large  variety  of  complex  curves  that 
cannot  be  modeled  by  either  the  standard  curves  or  the 
polynomials  of  the  conventional  design  methods.  Consequent- 
ly, B-splines  are  capable  of  solving  a  large  variety  of 
problems  and  produce  a  wider  range  of  solutions. 

The  integration  of  design  and  manufacture  has  also  been 
considered  for  the  case  of  plate  cams.  The  objectives  of 
considering  the  relevant  manufacturing  concerns  is  to  pro- 
duce a  design  which  requires  either  minimum  machining  time 
or  minimum  material  cost,  in  addition  to  satisfying  the 
specified  design  requirements. 
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It  is  a  tedious  procedure  to  manually  manipulate  B- 
splines  in  parametric  design  because  of  the  large  number  of 
design  parameters  involved.  In  order  to  overcome  this 
difficulty,  optimization  techniques  have  been  implemented  to 
automate  the  design  procedure. 

The  implementation  of  optimization  techniques  produced 
several  distinct  advantages.  These  techniques  successfully 
solved  the  test  problems  with  minimal  user  intervention, 
thus  reducing  the  amount  of  designer  time  and  the  skill 
level  required  to  produce  the  design.  They  also  allow  the 
designer  to  produce  the  best  design  with  respect  to  the 
specified  design  objective  and  performance  requirements. 

Optimization  techniques  also  allow  the  integration  of 
design  and  manufacture.  Integration  is  achieved  by  select- 
ing the  objective  function  based  on  the  manufacturing  con- 
cerns and  the  constraint  functions  based  on  the  performance 
requirements  that  the  design  must  satisfy.  The  objective 
function  is  either  the  machining  time  required  for  producing 
the  cam  on  an  NC  machine  or  the  material  cost. 

In  the  case  of  machining  time  minimization,  a  formula- 
tion to  estimate  the  cutting  time  was  derived  based  on  a 
variable  step  size  cutter  path  generation  algorithm.  In 
this  case,  the  cutting  time  was  used  as  the  cost  function. 
The  minimization  of  material  cost,  assuming  constant  depth, 
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is  equivalent  to  minimization  of  the  cross-sectional  area  of 
the  cam.  In  this  case,  the  cross-sectional  area  is  used  as 
the  cost  function  in  the  optimization  process. 

Several  constraints  were  considered  to  ensure  that  the 
final  design  is  safe  and  acceptable.  The  design  require- 
ments considered  were  limits  on  pressure  angle  and  contact 
stress,  along  with  elimination  of  separation,  oscillation, 
and  undercutting.  Other  design  requirements  can  be  conven- 
iently included  as  desired  through  user-supplied  constraint 
functions  which  can  be  input  through  user-supplied  subrou- 
tines. 

The  set  of  design  variables  considered  included  the 
base  radius  of  the  cam,  the  offset  of  the  roller  follower, 
the  orientation  of  the  cam  blank  on  the  NC  machine  table, 
and  the  coordinates  of  the  control  points  of  the  B-splines. 

An  exterior  penalty  function  method  was  used  to  solve 
the  optimization  problems  of  both  the  minimization  of  ma- 
chining time  and  the  minimization  of  material  cost.  For  the 
minimization  of  machining  time,  which  is  a  discontinuous 
function,  a  directed  grid  search  method  is  used  to  solve  the 
unconstrained  minimization  sub-problem.  For  the  minimiza- 
tion of  material  cost,  where  the  cost  and  constraint  func- 
tions are  continuous,  the  Fletcher-Reeves  conjugate-gradient 
method  is  used  for  unconstrained  minimization. 
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Several  numerical  examples  were  solved  to  verify  the 
proposed  design  technique.  In  all  the  examples,  this  tech- 
nique was  consistently  effective  in  constraint  correction 
and  cost  reduction.  The  design  procedure  was  highly  auto- 
mated and  required  little  designer  intervention.  The  method 
was  also  very  robust  and  able  to  handle  a  wide  variety  of 
problems,  which  indicates  the  great  potential  of  B-splines 
and  optimization  techniques  as  design  tools  in  cam  design. 

The  work  that  has  been  done  in  this  thesis  offers 
several  possibilities  for  future  development.  More  work  is 
needed  to  extend  this  technique  to  other  types  of  cam  and 
follower  systems.  It  is  also  recommended  that  this  tech- 
nique be  expanded  to  account  for  design  problems  involving 
cams  with  multiple  rise  cycles  and/or  multiple  return 
cycles.  Future  efforts  could  also  be  directed  towards 
treating  the  manufacturing  aspects  in  greater  detail.  It  is 
recommended  that  parameters  such  as  feed,  cutting  speed,  and 
tool  diameter  also  be  considered  as  design  variables  in  the 
optimization.  An  exploration  of  different  optimization 
techniques  to  improve  the  computational  efficiency  of  the 
design  process  is  also  necessary  for  enhanced  performance. 
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ABSTRACT 


Conventional  methods  of  designing  cams  are  very  re- 
strictive and  usually  produce  designs  which  only  guarantee 
satisfaction  of  the  boundary  conditions.   The  standard 
curves  and  polynomials  used  in  conventional  design  tech- 
niques do  not  provide  control  over  the  interior  of  a  motion 
cycle;  thus  the  design  obtained  may  violate  other  design 
requirements  such  as  pressure  angle  and  stress  limits.   The 
primary  objective  of  this  thesis  is  to  develop  an  automated 
design  technique  for  plate  cams  which  provides  control  over 
the  entire  profile  of  the  cam.   A  second  objective  is  to 
integrate  the  design  and  manufacture  of  plate  cams  in  order 
to  produce  designs  which  not  only  satisfy  the  design  re- 
quirements, but  also  ensure  efficiency  in  manufacture.  To 
achieve  the  first  objective,  B-spline  curves  are  used  to 
model  the  lift  diagram,  and  thereby  the  cam  profile.   The 
shape  of  the  interior  of  a  B-spline  is  only  partly  deter- 
mined by  the  boundary  conditions.   Thus,  the  curve  can  be 
designed  to  first  fit  the  boundary  conditions;  the  shape  of 
the  interior  can  then  be  modified  to  satisfy  other  require- 
ments without  affecting  the  curve  at  its  end  points.   Gener- 
ally, manual  design  of  a  B-spline  is  very  difficult  because 
of  the  large  number  of  the  design  variables  involved.   To 
overcome  this  difficulty,  a  design  technique  based  on  non- 
linear optimization  is  developed  for  the  purpose  of  automat- 
ing the  design  process.  Optimization  also  allows  the  design- 
er to  produce  the  best  design  for  a  problem  with  respect  to 
the  specified  requirements.   In  addition,  optimization 
techniques  allow  the  integration  of  design  and  manufacture 
by  permitting  requirements  from  both  aspects  to  be  included 
in  the  optimization  process.   Two  classes  of  problems  are 
considered  in  this  thesis:  design  for  minimum  machining  time 
and  design  for  minimum  material  cost.   In  both  cases,  the 
optimization  problems  are  solved  by  an  exterior  penalty 
function  method.   For  the  minimization  of  machining  time, 
which  is  a  discontinuous  function,  a  directed  grid  search 
method  is  used  to  solve  the  unconstrained  minimization  sub- 
problem.   For  the  minimization  of  material  cost,  where  the 
cost  and  constraint  functions  are  continuous,  the  Fletcher- 
Reeves  conjugate-gradient  method  is  used  for  unconstrained 
minimization.   The  design  technique  developed  in  this  thesis 
was  implemented  in  a  computer  program  and  tested  on  several 
examples.   In  all  the  examples,  the  design  technique  pro- 
duced good  designs  very  effectively.  The  design  process  was 
also  found  to  be  highly  automated  and  required  very  little 
designer  intervention.   It  may  be  concluded  that  B-splines 
are  a  powerful  tool  for  plate  cam  design  and  that  the  pro- 
posed optimization  technique  is  a  reliable,  robust,  and 
computationally  viable  method  for  automated  design  of  B- 
spline  cams. 


